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Soybean  [Glycine  max  (I.)  Merr.]  yield  has  been  shown  to  respond  to 
changes  in  plant  population.  Results  from  numerous  plant  population  studies 
on  soybean  conclusively  show  that  various  crop  and  plant  parameters  which 
are  known  to  influence  yield,  are  themselves  subject  to  plant  population  effects. 
However,  the  mechanism  by  which  plant  population  exerts  its  influence  on  yield 
Is  not  well  understood. 

Field  experiments  were  conducted  in  Gainesville,  Florida  (29°38N) 
during  1987  and  1988  to  identify  crop  parameters  and  physiological  and 
ecological  relationships  which  mediate  the  effect  of  plant  population  on  yield. 
SOYGRO--a  process-oriented  crop  growth  simulator,  was  also  used  to 
investigate  the  effect  of  plant  population  on  seed  yield  and  to  test  published 
hypotheses  regarding  the  efficiency  with  which  the  soybean  crop  uses 
photosynthetically  active  radiation  (PAR)  to  produce  seed  yield. 
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The  soybean  cultivars  Kirby  [MG  VIII],  Centennial  [MG  VI]  and  Forrest 
[MG  V],  were  planted  in  1987  at  populations  of  11.1,  25.0  and  44.4  plants  nr^  in 
a  square  arrangement.  The  parameters  recorded  at  growth  stages  R5,  R6  and 
R8  included  total  dry  weight,  pod  dry  weight  and  number,  and  seed  dry  weight 
and  number.  Plant  population  did  not  have  any  effect  on  the  variables 
recorded;  however,  all  variables  were  significantly  affected  by  cultivar.  Seed 
size  and  total  dry  weight  m"2  at  R5  accounted  for  the  differences  observed  in 
yield  between  the  cultivars. 

In  1987,  Kirby  [MG  VIII]  was  planted  at  populations  of  2.0,  4.0,  6.3,  11.1, 
16.0,  25.0,  44.4  and  69.4  plants  m-2  in  a  square  arrangement.  The  variables 
most  sensitive  to  plant  population,  and  which  were  also  closely  associated  with 
yield,  were  seed  number  and  pod  number  at  R8  and  total  dry  weight  and  leaf 
dry  weight  at  R5. 

The  SOYGRO  model  which  was  calibrated  to  give  responses 
approximating  those  of  observed  mean  yield  and  canopy  traits  to  plant  density, 
confirmed  the  hypothesis  that  seed  yield  efficiency  (seed  yield  /  unit  PAR 
intercepted)  tends  to  increase  with  plant  density.  This  increase  can  be 
attributed  in  part  to  an  increasing  ratio  of  mobilized-N  to  total-N  in  seeds  at 
higher  plant  densities. 
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CHAPTER  1 
INTRODUCTION 


Soybean  [Glycine  max  (L.)  Merr.],  the  prodigious  legume  which  was  first 
domesticated  approximately  3,000  years  ago  by  the  Chinese  has  come  to  be 
known  as  the  "yellow  jewel"  and  the  "wonder  bean."  Soybean  has  a  tremendous 
range  of  uses.  It  is  used  in  the  production  of  such  varied  products  as  cardboard, 
paint,  fire  extinguisher  foam  and  glue.  However,  because  of  its  high  protein  and 
oil  content,  soybean  finds  its  most  important  use  either  directly  or  indirectly  as 
food  for  both  humans  and  animals.  Soybean  is  viewed  by  some  as  a  major 
weapon  against  world  hunger.  But  if  the  full  potential  of  soybean  as  a  weapon 
against  hunger  is  to  be  realized,  the  processes  which  determine  its  grain  yield 
need  to  be  better  understood. 

Soybean  is  grown  under  diverse  cultural  conditions  and  usually  at  plant 
populations  which,  through  "  trial  and  error,  "  have  been  found  to  be  desirable. 
The  study  of  the  effect  of  plant  population  or  plant  population  density  on  soybean 
grain  yield  is  of  interest  for  two  main  reasons.  First,  for  a  given  cropping  system, 
the  quantitative  relationship  between  yield  and  plant  population  density  is 
important  in  estimating  i)  at  what  density  maximum  or  optimum  yield  can  be 
achieved  and  ii)  in  economic  terms,  which  plant  population  will  provide  the  most 
profitable  yield. 


Secondly,  plant  population  studies  allow  examining  the  process  of  yield 
achievement  and  its  determinants.  The  determinants  of  yield  or  potential  yield 
correlates  (PYCs)  can  be  defined  as  any  morphological  or  physiological 
characteristic  of  a  crop,  which  may  influence  final  seed  yield  at  one  time  or 
another  during  the  development  of  the  crop.  An  understanding  of  the  relationship 
between  PYCs  and  yield  is  basic  to  any  appreciation  of  the  process  of  yield 
achievement.  The  elucidation  of  the  nature  and  the  extent  of  the  functional 
relationship  between  PYCs  and  yield  is  facilitated  by  plant  population  studies.  A 
sufficiently  wide  range  of  plant  populations  may  provide  a  range  in  magnitude  of 
both  PYCs  and  yield,  and  this  allows  comprehensive  regression  analyses  to  be 
conducted  on  the  relationship  of  yield  to  PYCs. 

Therefore  plant  population  studies  are  important  in  understanding  both 
the  physiological  basis  of  yield  and  crop  improvement.  By  establishing  the 
relationships  i)  of  both  PYCs  and  yield  with  plant  population  density,  ii)  between 
PYCs  themselves  and  iii)  between  PYCs  and  yield,  a  basis  and  direction  for 
physiological  inquiry  into  the  process  of  yield  achievement  can  be  attained. 
Once  the  major  determinants  or  PYCs  of  yield  have  been  established,  these 
factors  may  also  serve  as  effective  selection  criteria  in  crop  improvement 
programs  which  seek  to  achieve  higher  yields,  and  foster  socioeconomic 
development. 

It  has  been  pointed  out  by  Duncan  (1986)  that  there  has  yet  to  emerge  a 
theory  capable  of  explaining  soybean  yield  within  the  context  of  plant  population 
studies.  In  an  attempt  to  appreciate  the  dynamics  of  yield  within  the  context  of 
plant  population  studies,  the  objectives  of  the  investigation  conducted  here  are; 
first,  to  develop  an  analytical  framework  for  the  analysis  of  yield-plant  population 


density  data;  secondly,  to  ascertain  which  biological  and  ecological  variables 
and  relationships  of  the  soybean  crop  mediate  the  influence  of  plant  population 
density  on  final  grain  yield;  and  thirdly,  to  elucidate  the  basis  for  the  response 
surface  which  characterizes  yield-plant  population  density  relationships.  This 
final  objective  which  focuses  on  a  relatively  unexplored  area  (the  response 
surface)  of  plant  population  studies,  will  be  carried  out  using  SOYGRO  -  a 
process  oriented  crop  growth  simulation  model  for  soybean. 


CHAPTER  2 
LITERATURE  REVIEW 


The  research  literature  on  yield-plant  population  density  studies  will  be 
broken  down  into  three  sections.  First,  the  effect  of  plant  population  on  the 
soybean  crop  will  be  discussed.  The  next  two  sections  come  under  the  umbrella 
of  explaining  yield,  and  consist  of  discussions  on  the  observation  that  yield  per 
unit  area  increases  and  yield  per  plant  decreases  with  increasing  plant 
population  density,  and  discussions  on  the  asymptotic  or  parabolic  shape  of 
yield-plant  population  density  responses.  It  is  within  the  context  of  the  latter 
section  that  SOYGRO-a  process  oriented  plant  growth  simulator  (Jones  et  al., 
1989)  will  be  discussed.  Of  specific  interest  with  SOYGRO  are  estimations  of  the 
impact  of  plant  population  density  on  the  efficiency  with  which  plants  use  light 
inputs  for  seed  production.  Lastly,  the  identification  of  the  reproductive  growth 
stages  of  soybean  and  the  square  planting  pattern  which  are  relevant  to 
experimental  procedures  will  be  discussed. 

The  Effect  of  Plant  Population  on  Soybean 

In  plant  population  or  yield-plant  population  density  studies,  the  main 
ecological  force  which  influences  crop  performance  is  interplant  competition. 
With  increasing  density,  plants  tend  to  compete  more  intensely  for  available 


resources.  Plant  competition  is  subject  to  a  number  of  differing  interpretations; 
however,  here  it  will  be  defined  according  to  Grime  (1973)"the  tendency  for 
neighboring  plants  to  utilize  the  same  quantum  of  light,  ion  of  a  mineral,  molecule 
of  water  or  volume  of  space. 

Competitive  interactions  between  organisms  are  considered  to  be  among 
the  most  important  ecological  phenomena  (Snell  and  Burch,  1975).  Competition 
between  organisms  occupies  a  place  of  prominence  in  evolutionary  theory, 
where  "survival  of  the  fittest"  determines  the  procreation  possibilities  of  species. 
At  a  more  immediate  level,  competition  holds  various  consequences  for  a 
developing  plant.  In  reacting  or  adapting  to  changes  in  the  pressures  of 
competition,  plants  may  become  modified  both  physiologically  and 
morphologically.  When  this  occurs,  the  plant  is  said  to  exhibit  plastic  responses 
(Hutchings  and  Budd,  1981).  Plastic  responses  may  affect  the  productive 
performance  of  plants  and  crops.  It  has  been  amply  demonstrated  by 
researchers,  that  in  soybean,  both  yield  and  PYCs  are  subject  to 
density-mediated  plastic  responses.  The  impact  of  plant  population  density  on  i) 
vegetative  and  reproductive  structures,  ii)  lodging  and  iii)  nitrogen  (N)  and 
phosphorus  (P)  accumulation  will  be  discussed. 

Stems  and  Branches 

The  height  and  thickness  of  soybean  stems  have  been  shown  to  be 
affected  by  plant  population  density.  It  has  been  reported  that  plant  height 
increases  with  plant  population  density  (Ramseur  et  al.,  1985;  Doss  and  Thurlow, 
1974;  Wilcox,  1974;  Johnson  and  Harris,  1967).  This  increase  in  plant  height  has 


been  attributed  to  an  increase  in  the  length  of  the  internodes,  and  not  to 
increases  in  the  number  of  nodes  (Dominguez  and  Hume,  1978;  Basnet  et  al., 
1974;  Hinson  and  Hansen,  1962).  In  some  cases,  plant  height  has  been 
unaffected  by  density  (Hoggard  et  al.,  1978;  Lueschen  and  Hicks,  1977;  Probst, 
1945).  In  one  atypical  case,  Dominguez  and  Hume  (1978)  reported  that  with  the 
cultivar  'Vansoy',  increased  density  resulted  in  shorter  plants.  The  diameter  of 
the  main  stem  in  soybean  decreases  with  increased  density  (Miura  and  Gemma, 
1986;  Ramseur  et  al.,  1985;  Wright  et  al.,  1984;  Pontes  and  Ohirogge,  1972). 
Oba  et  al.  (1961)  also  found  that  stem  weight  per  plant  decreases  with  density. 

The  number  of  branches  per  plant  decreases  with  density  (Basnet  et  al., 
1974;  Pontes  and  Ohirogge,  1972;  Hinson  and  Hansen,  1962;  Lehman  and 
Lambert,  1960).  However,  in  their  work  Chaudhry  and  Cheema  (1985)  found 
that  the  effect  of  density  on  the  number  of  branches  per  plant  was  not  significant. 

Leaves  and  Leaf  Area  Index 

Leaf  weight  per  plant  decreases  (Oba  et  al.,  1961)  while  the  leaf  area 
index  (LAI)  increases  with  plant  population  density  (Parvez  et  al.,  1989;  Weber  et 
al.,  1966).  Of  relevance  to  population  studies  is  the  fact  that  the  degree  of 
shading  or  exposure  to  light  affects  leaf  anatomy.  At  lower  light  intensities  leaves 
develop  fewer  layers  of  palisade  cells  (Xu  and  Miao,  1988),  while  higher  light 
intensities  have  the  opposite  effect,  resulting  in  thicker  leaves  (Chabot  et  al., 
1979;  Pails  et  al.,  1982).  At  higher  population  densities  where  the  shading  of 
lower  leaves  occurs,  it  can  be  expected  that  there  will  be  greater  differences 


between  the  specific  leaf  weight  of  leaves  on  the  periphery  of  the  canopy  and 
those  deeper  in  the  canopy  itself. 

Flower  Production  and  Development 


Increased  population  results  in  a  reduction  in  the  proportion  of  flowers 
which  result  in  mature  pods  (Buttery,  1969).  Dominguez  and  Hume  (1978)  also 
found  that  increased  density  increased  the  percentage  of  flowers  aborted.  They 
further  observed  that  the  decrease  in  the  number  of  flowers  per  plant  observed  at 
high  densities,  was  due  to  fewer  flowers  produced  per  node.  This  observation 
has  also  been  made  by  Lehman  and  Lambert  (1960).  Fewer  flowers  at  higher 
densities  has  also  been  attributed  to  fewer  nodes  per  plant  (Dominguez  and 
Hume,  1978). 

Nodes.  Pods  and  Seeds 

The  number  of  nodes  per  plant  decreases  with  plant  population  density 
(Dominguez  and  Hume,  1978;  Basnet  et  al.,  1974).  However,  on  a  per  unit  area 
basis,  the  number  of  nodes  increases  with  density  (Parvez  et  al.,  1989;  Hruska 
and  Labounek,  1985).  Consistent  with  these  observations  are  those  that  the 
number  of  pods  per  plant  decreased  (Olsen,  1986;  Nakagawa  et  al.,  1986; 
Hoggard  et  al.,  1978;  Lueschen  and  Hicks,  1977;  Basnet  et  al.,  1974;  Lehman 
and  Lambert,  1960),  while  the  number  of  pods  per  unit  area  increased  (Udoguchi 
and  McCloud,  1986)  with  plant  population  density.  In  their  work,  Chaudhry  and 
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Cheema  (1985)  reported  that  the  number  of  pods  per  plant  was  not  affected  by 
plant  population  density.  Pod  set  is  also  greater  at  lower  populations  (Oba  et  al., 
1961;  Weber  etal.,  1966). 

Seed  weight  per  unit  area  Increases  with  plant  population  density  (Parvez 
et  al.,  1989;  Udoguchi  and  McCloud.  1986;  Parker  et  al.,  1981;  Spilde  et  a!., 
1980;  Doss  and  Thurlow  1974).  However,  in  their  work,  Hoggard  et  al.  (1978) 
found  that  seed  weight  per  unit  area  decreased  with  increasing  plant  population 
density.  Basnet  et  al.  (1974)  and  Caviness  (1966)  reported  that  seed  weight  per 
unit  area  was  unaffected  by  plant  population  density. 

Johnson  and  Harris  (1967)  in  their  work,  found  that  of  the  several  cultivars 
investigated,  weight  per  seed  decreased  with  plant  population  density  only  with 
'Lee'.  This  decrease  has  also  been  reported  by  Parks  and  Manning  (1980), 
Wilcox  (1974),  and  Wright  et  al.  (1987).  In  thinning  experiments,  Weil  and 
Ohirogge  (1976)  also  reported  that  weight  per  seed  increased  by  11.5%  with 
decreased  plant  population  density.  Nedic  et  al.  (1987)  also  found  that  the 
number  of  filled  seeds  per  plant  decreased  with  density.  The  reduction  in  seed 
size  with  increasing  plant  population  density  may  be  the  result  of  increased  pod 
production  per  unit  area.  Other  researchers  reported  no  relationship  between 
weight  per  seed  and  plant  population  density  (Dominguez  and  Hume,  1978; 
Pontes  and  Ohirogge,  1972). 

The  number  of  seeds  per  plant  decreases  with  density  (Lueschen  and 
Hicks,  1977;  Pontes  and  Ohirogge,  1972).  Blumenthal  et  al.  (1988)  observed 
that  the  number  of  seeded  pods  per  branch  decreased  with  population.  The 
number  of  seeds  per  pod  decreases  with  plant  population  density  (Blumenthal  et 
al.,  1988;  Herbert  and  Litchfield,  1982;  Pontes  and  Ohirogge,  1972;  Lehman  and 


Lambert,  1960).  However,  in  some  cases  it  has  been  shown  that  seeds  per  pod 
is  not  affected  by  plant  population  density  (Chaudhry  and  Cheema,  1985; 
Dominguez  and  Hume,  1978). 

The  distribution  of  pods  and  hence  seeds  within  the  crop  and  on  plants  is 
also  affected  by  plant  population  density.  The  height  of  the  first  pod  increases 
with  density  (Nedic  et  al.,  1987;  Olsen,  1986;  Tsuchiya  et  al.,  1986;  Lueschen  and 
Hicks,  1977;  Basnet  et  a!.,  1974).  This  fact  is  of  importance  in  mechanical 
harvesting,  where  it  is  advantageous  to  have  the  pods  located  away  from  the 
ground.  Parks  and  Manning  (1980)  found  that  the  ratio  of  the  number  of  seeds  on 
the  main  stem  to  those  on  the  branches  increased  with  plant  population  density. 
Ramseur  et  al.  (1984)  also  found  that  with  decreased  plant  population  density, 
the  seed  number  and  yield  on  branches  increased  while  pod  number,  seed 
number  and  yield  on  stems  decreased.  Dominguez  and  Hume  (1978)  in  their 
work,  observed  that  for  the  determinate  cultivar  'Fiskeby  V,  an  increase  in  plant 
population  density  resulted  in  a  shift  in  the  majority  of  yield  from  the  upper  third  to 
the  middle  of  the  plant.  Ramseur  et  al.  (1984)  also  found  that  increased  plant 
population  density  reduced  the  contribution  of  the  lower  nodes  to  yield. 
According  to  Weil  and  Ohirogge  (1976)  thinning  resulted  in  weight  per  seed 
being  23.9%  greater  at  the  top  of  the  canopy  than  in  unthinned  stands.  Increased 
density  resulted  in  a  reduction  of  the  bottom  portion  of  plants  in  their  contribution 
to  total  seed  number  (Dominguez  and  Hume,  1978). 


10 
Biomass 

It  has  been  reported  that  biomass  per  unit  area  increases  with  plant 
population  density  (Parvez  et  a!.,  1989;  Blumenthal  et  al.,  1988;  Udoguchi  and 
McCloud,  1986).  Biomass  per  plant  decreases  with  plant  population  density 
(Holler  and  Abrahamson,  1977;  Wilcox,  1974).  It  should  be  noted  that  biomass 
per  unit  area  tends  to  increase  with  plant  population  density  at  all  stages  of  plant 
development.  However,  biomass  per  plant  will  decrease  with  plant  population 
density  only  when  there  is  appreciable  interplant  competition. 

Lodging 

Lodging  increases  with  plant  population  density  (Tsuchiya  et  al.,  1986; 
Pontes  and  Ohirogge,  1972;  Probst,  1945).  Lueschen  and  Hicks  (1977)  found  a 
significant  linear  relationship  between  plant  population  and  lodging.  Johnson 
and  Harris  (1967)  in  studying  a  number  of  soybean  cultivars  reported  that 
cultivars  with  larger  stem  diameters,  such  as  'Hardee',  tended  not  to  lodge  in 
response  to  increased  plant  population  density.  Increased  lodging  was  also 
associated  with  taller  plants  (Cooper,  1971).  Increased  density  increased  plant 
height  and  decreased  stem  diameter  (Wright  et  al.,  1984). 

Nitrogen  and  Phosphorus  Accumulation 

According  to  Abrahamson  and  Caswell  (1982),  biomass  allocation  in 
plants  has  been  correlated  to  various  degrees  to  the  allocation  of  minerals  and 
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energy.  Hanway  and  Weber(1971)  from  their  work,  also  showed  that  N,  P  and  K 
accumulation  mirrored  biomass  allocation.  Nelson  and  Weaver  (1980)  reported 
that  the  lowest  planting  density  exhibited  the  highest  rate  of  acetylene  reducing 
activity  per  plant.  This  finding  has  been  attributed  to  the  greater  nodule  mass  on 
plants  at  the  lowest  plant  population  density,  since  the  specific  activity  of  nodules 
was  only  slightly  affected  by  plant  population  density.  Buttery  (1969)  reported  that 
N  concentration  decreased  with  density,  and  he  attributed  this  to  increased 
competition  for  N  between  plants  at  higher  plant  population  densities. 

The  nutrient  status  of  plants  is  directly  related  to  crop  productivity  in  two 
ways.  First,  the  photosynthetic  capacity  of  plant  tissue  depends  upon  the 
concentration  of  nutrients,  especially  in  the  leaves.  Secondly,  the  status  of  N 
and  P  also  determines  the  amount  of  nutrient  which  may  be  available  for  seed 
production  via  the  process  of  nutrient  remobilization. 

The  Structural  and  Yield  Implications  of  Density-Mediated  Plastic  Responses 

Plastic  responses,  which  are  embodied  in  morphological  and 
physiological  changes,  hold  both  structural  and  yield  implications  for  soybean. 
Structurally,  plants  are  exposed  to  various  physical  stresses.  Plants  which 
undergo  density-mediated  plastic  responses  must  compensate  for  changing 
stresses  to  avoid  the  mechanical  failure  of  their  axes  (that  is,  breakage  and 
lodging).  Accordingly,  plants  preferentially  form  and  lay  down  specialized  tissues 
of  varying  material  properties,  in  different  regions  of  the  plant  (Niklas  and 
O'Rouke,  1982).  Pianka  (1972)  suggested  that  in  organisms  with  "plastic 
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abilities,"  it  can  be  expected  that  the  manner  and  pattern  of  resource  allocation 
will  be  influenced  by  different  competitive  environments. 

The  process  of  adjusting  to  the  forces  of  interplant  competition,  therefore, 
influences  canopy  arrangement,  light  interception,  stem  thickness,  pod  number 
and  distribution,  mineral  concentrations  and  contents,  the  amount  of  nutrients 
which  can  be  remobilized,  the  number  of  seeds  and  so  forth.  The  entire  process 
whereby  light  energy  is  converted  into  chemical  energy  and  is  allocated  and 
stored,  is  affected  by  plant  population  density. 

Explaining  Soybean  Yield  in  Plant  Population  Studies 

In  attempting  to  explain  soybean  yield  within  the  context  of  plant 
population  studies,  three  areas  of  overlapping  research  can  be  identified:  i)  the 
correlation  or  regression  of  PYCs  with  yield,  ii)  the  exploration  of  the 
physiological  basis  for  the  observation  that  yield  per  unit  area  increases,  while 
yield  per  plant  decreases  with  increasing  plant  population  density,  and  iii)  the 
examination  of  the  physiological  basis  of  asymptotic  and  parabolic  yield-plant 
population  density  responses. 

Regression  Analyses  and  Correlation  of  PYCs  and  Yield 

The  correlation  or  association  of  various  plant  and  crop  characteristics 
with  seed  yield  is  standard  practice  in  yield  improvement  programs  and  in  yield 
analysis.  Outside  the  context  of  plant  population  density  studies  (using  different 
soybean  genotypes  and  productive  environments),  the  following  characteristics 
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are  correlated  with  yield  per  plant  In  soybean:  stem  weight,  stem  thickness  and 
branch  number  (Taguchi  et  al.,  1958),  pod  number  and  the  number  of  fertile 
nodes  (Board,  1985).  Yield  per  unit  area  has  been  shown  to  be  highly 
associated  with  LAI  and  number  of  seeds  per  unit  area,  but  only  moderately 
correlated  with  dry  matter  production  (Basuchaudhuri,  1987).  In  their  work,  Kaw 
and  Menon  (1972)  found  that  pod  number  and  seed  number  per  unit  area  were 
highly  correlated  with  yield  per  unit  area.  Characteristics  showing  the  lowest 
correlation  with  yield  per  plant  include  seeds  per  pod  and  weight  per  seed 
(Board,  1985;  Tandojam,1986;  Johnson  et  al.,  1955).  Hartwig  and  Edwards 
(1970)  in  their  studies  on  morphological  characters  which  affect  yield  per  unit 
area,  reported  that  seed  number  per  pod  did  not  affect  yield.  However,  Ma 
(1946)  (cited  by  Anand  and  Torrie  (1963))  found  that  seed  number  per  pod  was 
correlated  with  yield,  and  Weber  and  Moorthy  (1952)  reported  that  there  was  a 
high  positive  correlation  between  weight  per  seed  and  yield  per  unit  area. 

Within  population  studies  the  following  characteristics  are  highly 
correlated  with  yield  per  plant:  stem  diameter  and  pod  number  per  plant  (Pontes 
and  Ohirogge,  1972).  The  following  were  correlated  with  yield  per  unit  area:  pod 
number  and  total  dry  matter  (Parvez  et  al.,  1989;  Udoguchi  and  McCloud,  1986). 
Leaf  area  index  was  not  correlated  with  yield  per  unit  area  (Ramseur  et  al.,  1985; 
Costa  etal.,  1980). 


Plant  Population  Density  Effect  on  Productivitv  on  a  per  Plant  and  per  Unit  Area 
Basis 


Light  interception  and  canopy  structure  form  the  main  basis  of  discussions 
which  seek  to  explain  the  impact  of  plant  population  density  on  yield  per  plant 
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and  yield  per  unit  area.  To  a  lesser  extent,  the  significance  of  N  and  P  status  has 
also  been  addressed.  Writing  in  1966,  Weber  et  al.  argued  that  in  order  to 
understand  how  plants  respond  to  changes  in  density,  it  is  necessary  to 
investigate  the  'spatial  relationships'  in  canopies,  solar  energy  interception  and 
leaf  area  accumulation.  The  amount  of  light  intercepted  by  a  crop  is  determined 
by  the  arrangement  of  the  canopy  (branching,  leaf  area,  etc.) 

Lueschen  and  Hicks  (1977)  reported  that  soybean  plants  are  capable  of 
compensating  for  low  plant  population  densities  by  producing  more  branches 
and  pods  per  plant.  Acock  and  Acock  (1987)  found  that  the  growth  and 
development  of  branches  on  a  plant  depended  on  the  amount  of  light  received. 
Shading  reduced  the  photosynthetic  rate  of  plants  at  high  densities,  resulting  in 
carbon  shortages  (per  plant)  and  a  tendency  to  produce  fewer  branches.  These 
findings  are  consistent  with  those  of  Charles-Edwards  (1984),  Charles-Edwards 
and  Beech  (1984),  Bunce  (1988)  and  Johnston  et  al.  (1969). 

Yield  per  unit  area  has  been  positively  correlated  with  vegetative  dry 
matter  production  per  unit  area  (Egli  et  al.,  1987).  It  has  also  been  established 
that  dry  matter  production  is  a  function  of  solar  radiation  intercepted  (Shibles  and 
Weber,  1966;  Brougham,  1956;  Davidson,  1958;  Watson,  1958).  Brougham 
(1956)  found  that  both  the  rate  of  dry  matter  production  per  unit  area  and  the 
percent  solar  radiation  intercepted  increased  with  LAI.  Johnston  and  Pendleton 
(1968)  argued  that  due  to  its  large  leaves  and  phyllotaxy,  soybeans  tend  to  have 
a  closed  canopy.  Soybeans  intercepted  as  much  as  90%  of  the  incoming  solar 
radiation  near  the  top  and  periphery  of  the  canopy  (Sakamoto  and  Shaw,  1967). 
Thus  there  is  a  tendency  for  lower  leaves  to  be  shaded.  If  increased  plant 
population  density  is  accompanied  by  increased  shading,  it  can  be  expected  that 
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the  light  intercepted  per  plant  will  decrease  which  will  subsequently  result  in  a 
reduction  in  yield  per  plant. 

Leaf  area  index  increases  with  plant  population  density  (Parvez  et  al., 
1989;  Weber  et  al.,  1966;  Costa  et  al.,  1980).  The  extent  of  photosynthetic  activity 
not  only  depends  on  LAI  and  the  amount  of  light  intercepted,  but  also  upon  the 
photosynthetic  capacity  of  the  leaves  intercepting  the  light.  Bunce  (1988)  found 
that  photosynthetic  capacity  was  lower  at  higher  densities,  even  though  LAIs  are 
higher.  Reduced  photosynthetic  capacity  resulted  from  the  shading  of  lower 
leaves,  and  the  longer  period  of  time  that  young  developing  leaves  spent  in  the 
shade  at  high  stand  densities.  These  findings  corroborate  the  work  of 
Nichiporovich  et  al.  (1969)  and  the  suggestion  by  Johnson  et  al.  (1969)  that  due 
to  both  inter-  and  intra-plant  competition  for  light,  middle  and  bottom  soybean 
leaves  do  not  reach  their  full  photosynthetic  potential. 

Differences  in  light  intensity  also  results  in  photomorphogenic  changes  in 
leaves.  Phytochrome  provides  the  sensory  pigment  which  enables  plant  cells  to 
react  to  light  quality  and  intensity.  The  light-regulated  ratio  of 
phytochrome-far-red  to  total  phytochrome  serves  to  influence  protein  synthesis 
and  enzyme  activity  (Song  1984;  Schopfer  1984).  Photomorphogenic  effects 
have  implications  for  photosynthetic  activity.  Increased  irradiance  results  in 
thicker  leaves,  and  this  has  been  attributed  to  increased  thickness  of  the  palisade 
mesophyil  layer  (Reginer  et  al.,  1988;  Chabot  et  al.,  1979).  Therefore, 
photosynthetic  tissue  per  unit  leaf  area  is  increased.  Ford  et  al.  (1983)  found 
that  carbon  dioxide  exchange  rate  was  moderately  and  positively  correlated  with 
specific  leaf  weight.  Wells  et  al.  (1986a)  also  found  that  differences  in  canopy 
photosynthesis  between  soybean  cultivars  was  associated  with  differences  in 
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specific  leaf  weight.  However,  by  contrast,  Secor  et  al.  (1982)  found  that  specific 
leaf  weight  was  able  to  explain  only  4%  of  the  variation  in  total  photosynthesis 
among  110  lines  of  soybean.  Nelson  and  Schweitzer  (1988)  found  that  specific 
leaf  weight  in  soybean  is  an  effective  selection  trait,  as  it  tends  not  to  be  greatly 
influenced  by  the  environment. 

Wells  et  al.  (1986b)  reported  that  the  quantity  of  RUBPcase  per  unit  leaf 
area  was  positively  correlated  with  specific  leaf  weight.  Reginer  et  al.  (1988) 
reported  that  when  soybean  plants  were  grown  at  low  versus  high  irradiance,  the 
total  soluble  protein  per  unit  leaf  area  decreased  by  as  much  as  40%,  RuBPcase 
per  unit  area  decreased  by  31%,  and  chlorophyll  content  per  unit  leaf  area 
decreased  by  15%.  These  factors  contribute  to  reduced  photosynthetic  rates  at 
lower  irradiance.  However,  in  soybean,  shaded  leaves  are  not  'parasitic,' 
because  their  respiration  rates  decrease  along  with  reductions  in  photosynthetic 
capacity  (Shibles  and  Weber  1965;  Duncan  et  al.,  1967). 

The  status  of  N  in  both  leaf  tissue,  and  the  entire  plant,  is  important  to  the 
processes  of  photosynthesis  and  remobilization,  both  of  which  are  related  to 
seed  yield.  The  rate  of  carbon  dioxide  assimilation,  chlorophyll  content  and 
RuBPcase  activity  are  proportional  to  leaf  N  (Evans,  1983).  It  has  been  shown  by 
Boote  et  al.  (1978)  that  a  decline  in  photosynthesis  was  correlated  with 
decreases  in  leaf  N.  Lugg  and  Sinclair  (1981)  reported  that  in  lower  leaflets 
and  below  a  critical  protein  content  per  unit  leaf  area  (13  to  16  g  m'^)  in 
uppermost  leaflets,  there  existed  a  linear  relationship  between  estimated  protein 
content  per  unit  area  and  photosynthesis.  Buttery  (1969)  found  that  N 
concentration  decreased  with  increasing  plant  population  density  and  attributed 
this  to  an  inadequate  supply  of  N.  Soluble  protein  and  N  content  in  the  leaf  reach 
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a  maxima  during  the  early  reproductive  phase  (Boon-Long  et  al.,  1983, 
Thibodeau  and  Jaworski,  1975).  Pod  number  and  seed  yield  are  correlated  with 
N  accumulation  in  leaves  in  the  early  reproductive  phase  (Zheng  et  al.,  1987). 

Phosphorus  status  also  affects  the  productive  capacity  of  soybean 
canopies.  Fredeen  et  al.  (1989)  found  that  low  P  decreased  soybean  growth 
primarily  through  effects  on  the  rate  of  expansion  of  leaf  surfaces.  Photosynthetic 
rates  were  also  reduced  (Fredeen  et  al.,  1989).  Sawada  et  al.  (1983)  in  their 
P-stress  experiments  reported  that  decreases  in  photosynthetic  activity  can  be 
partially  attributed  to  P  deficiency  (lack  of  phosphate  compounds  in  the 
photosynthetic  machinery).  Phosphorus  deficiency  reduced  photosynthetic 
capacity  (Foyer  and  Cooper,  1986).  Photosynthetic  CO2  fixation  diminished  with 
reduced  P  concentration  in  leaves  (Terry  and  Ulrich,  1973).  Keogh  et  al.  (1972) 
found  that  P  levels  in  10  soybean  cultivars  from  three  maturity  groups  ranged 
from  3.4  to  3.7  g  kg'"' ,  with  low  standard  deviations.  Nutrient  concentration  did 
not  vary  much  except  Ca  which  was  higher  in  higher-yielding  cultivars. 

Miller  et  al.  (1961)  reported  that  in  the  canopy  of  soybean,  upper  leaves 
tend  to  have  a  higher  P  concentration  than  lower  leaves.  They  also  found  a 
positive  relationship  between  P  concentration  (in  upper  leaves)  and  yield  when 
percent  K  was  greater  than  1.2  g  kg'"".  Saturating  irradiance  for  photosynthesis 
in  P-deficient  plants  was  found  to  be  30%  less  than  in  the  control,  however,  P 
deficiency  did  not  affect  photosynthetic  activity  at  low  irradiances  (Sawada  et  al., 
1983).  Thus  at  higher  plant  population  densities  (if  there  are  no  deficiencies  due 
to  interplant  competition)  low  P  would  affect  leaves  atop  the  canopy. 

Photosynthetic  capacity  and  seed  production  are  also  affected  by  the 
process  of  nutrient  remobilization.  Nitrogen  and  P  are  accumulated  in  the  leaves. 
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petioles  and  stems  until  the  pod  production  stage,  after  which  they  are  rapidly  lost 
from  these  parts  and  accumulated  in  the  developing  seeds  (Hanway  and  Weber, 
1971).  In  upper  leaves,  N  and  P  concentrations  were  50.5  and  3.7  g  kg"'' 
respectively,  while  in  the  bottom  leaves  (older  and  undergoing  translocation)  N 
and  P  concentrations  were  23.0  and  3.0  g  kg"''  respectively  (Terman,  1977). 

The  extent  to  which  developing  fruits  draw  upon  nutrient  reserves  of 
vegetative  tissue  has  led  to  the  suggestion  that  nutrient  remobilization  leads  to 
senescence.  Depodding  prevented  senescence  and  the  gradual  reinstatement 
of  senescence  as  an  increasing  number  of  pods  is  allowed  to  develop  on  each 
plant  (Lindoo  and  Nooden,  1977).  However,  Nooden  et  al.  (1978)  indicate  that 
the  development  of  seeds  is  separable  from  the  senescence  response  in 
soybeans  and  that  seeds  may  function  as  more  than  sinks.  They  found  that  sink 
size  does  not  parallel  foliar  senescence;  the  senescence  response  is  saturated  at 
a  level  far  below  the  maximum  level  of  dry  weight  and  nitrogen  accumulation  in 
the  seeds.  Although  drain  and  diversion  may  be  involved  in  the  monocarpic 
senescence  of  soybeans,  it  seems  unlikely  that  developing  seeds  exert  such  a 
remarkable  correlative  influence  simply  by  functioning  as  a  sink  (Nooden  et  al., 
1978). 

Carbon  assimilated  early  in  the  life  cycle  has  virtually  no  direct  relevance 
to  fruit  nutrition.  Therefore,  carbon  fixed  during  reproductive  development  must 
be  primarily  relied  upon  to  furnish  the  seed's  requirement  for  this  element.  The 
extremely  low  rate  of  transfer  of  early-fixed  carbon  to  seeds  is  resolved  as  being 
partly  due  to  the  fact  that  much  of  this  carbon  has  been  dissipated  in  respiration 
before  flowering  has  commenced,  and  partly  to  the  fact  that  carbon  which 
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survives  until  fruiting  is  bound  into  materials  from  which  it  cannot  be  readily 
retrieved  during  the  senescence  of  vegetative  parts  (Pate  and  Flinn,  1973). 

On  the  other  hand,  Pate  and  Flinn  found  that  N  assimilated  early  in  the 
life  cycle  is  released  for  seed  development  with  great  efficiency.  Therefore,  it 
should  be  theoretically  possible  to  increase  rates  of  transfer  to  seeds  by 
arranging  for  larger  reserves  of  N  to  have  accumulated  by  the  time  when 
flowering  commences.  Mobilization  is  a  gradual  process  which  gathers 
momentum  as  fruiting  proceeds.  Hocking  and  Pate  (1977)  in  studying  legumes 
reported  that  N,  P  and  K  are  highly  mobile  and  60  to  90%  are  usually  retrieved. 
Approximately  half  of  the  N  and  P  in  mature  seeds  appears  to  have  been  in  other 
plant  parts  prior  to  pod  development  (Hanway  and  Weber,1971).  Therefore, 
since  nutrient  accumulation  mirrors  biomass  accumulation,  at  higher  plant 
population  densities  more  substrate  will  be  available  per  unit  area  for 
remobilization   making  higher  seed  yields  possible. 

The  Yield-Plant  Population  Density  Response 

Explanations  of  the  dynamics  of  yield  in  population  studies,  go  beyond 
justifying  that  yield  per  area  tends  to  increase  while  yield  per  plant  tends  to 
decrease  with  plant  population  density.  The  asymptotic  or  parabolic  shape  of  the 
yield-plant  population  density  response  needs  to  be  investigated.  The  main 
areas  of  interest  in  explaining  the  shape  of  yield-plant  population  density 
reponses  are  light  interception,  the  efficiency  with  which  the  crop  and  plant  use 
production  inputs,  biomass  accumulation,  lodging  and  barreness.  In  this 
subsection,  first,  the  fundamental  and  obvious  characteristics  of  the  yield-plant 
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population  density  response  will  be  discussed;  secondly,  Duncan's  (1986) 
theory  on  yield-plant  population  density  response  in  soybean  will  be  critically 
reviewed;  thirdly,  the  relevance  of  SOYGRO  to  investigations  of  yield-plant 
population  density  responses  will  be  outlined;  and  fourthly,  the  possible  basis  for 
changes  in  the  efficiency  with  which  plants  use  inputs  for  seed  production  and  its 
relevance  to  yield-plant  population  density  responses  will  be  discussed. 

The  fundamentals  of  yield-plant  population  density  responses 

In  Fig.  2-1  (a)  and  Fig.  2-1  (b)  the  typical  yield  -  plant  population  density 
responses  are  shown.  Using  plant  population  density  as  a  basis  for  categorizing 
sections  of  these  responses,  the  following  phases  are  identified:  between  plant 
population  densities  DO  and  D1 ,  yield  per  plant  is  at  a  maximum  and  is  constant, 
and  yield  per  area  assumes  a  strict  linear  relationship  with  plant  population 
density.  At  density  D1,  interplant  competition  commences  and  yield  per  plant 
begins  to  decrease.  After  plant  population  density  D1 ,  yield  per  area  begins  to 
deviate  from  linearity  and  the  slope  of  the  response  begins  to  decrease. 
Between  plant  population  densities  D1  and  D4,  although  yield  per  plant 
decreases,  yield  per  area  increases  because  plant  population  density  increases 
at  a  rate  greater  in  magnituide  than  the  rate  of  reduction  in  yield  per  plant.  At 
plant  population  density  D2  the  yield  per  plant  response  assumes  its  steepest 
negative  slope,  and  between  plant  population  densities  D2  and  D3  (which  occurs 
at  an  early  stage  of  interplant  competition)  interplant  competition  has  its  greatest 
impact  on  the  reduction  of  yield  per  plant.  This  impact  is  given  by  the  magnitude 
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Fig.  2-1.  Typical  responses  for  yield  per  unit  area  (a)  and  yield  per 
plant  (b)  to  plant  density. 
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of  the  negative  slope  of  the  relationship  between  yield  per  plant  and  plant 
population  density  between  D2  and  D3. 

At  plant  population  density  D3,  there  is  a  reduction  in  the  rate  of  decrease 
in  yield  per  plant  with  plant  population  density.  This  reduction  may  be  indicative 
of  plant  adaptation  to  interplant  competition,  and  possibly  implies  an  increased 
efficiency  with  which  the  plant  is  using  inputs  for  seed  production. 

Between  plant  population  densities  D4  and  D5  yield  per  area  assumes  a 
maximum  and  constant  value.  A  number  of  crop  and  plant  requirements  for 
achieving  maximum  yield  per  area  have  been  suggested.  It  has  been  suggested 
that  insolation  interception  must  approach  100%  early  in  reproductive  growth  to 
maximize  yield  per  area  (Shibles  and  Weber,  1965;  Hawkins,  1982),  specifically, 
prior  to  the  period  of  grain  production  (Shibles  and  Weber,  1966).  Duncan 
(1986)  argued  that  complete  interception  needed  to  occur  only  before  fruit  or 
seed  number  reached  a  maximum,  in  order  to  maximize  fruits  per  unit  area. 
Duncan's  view  may  be  consistent  with  the  findings  of  Egli  et  al.  (1985)  that  fruit 
number  increased  even  after  growth  stage  R5.  However,  Johnson  et  al.  (1982) 
have  suggested  that  complete  light  interception  is  required  at  the  earlier  growth 
stage  of  R1 ,  if  yield  is  to  be  maximized. 

Vegetative  mass  has  also  been  correlated  to  seed  number.  This 
observation  reveals  the  possible  importance  of  vegetative  dry  weight  in 
contributing  to  differences  in  yield  where  light  interception  is  constant,  as  pointed 
out  by  Duncan  (1986).  Egli  et  al.  (1987)  reported  that  increasing  the  vegetative 
dry  weight  to  a  limit  of  500  to  600  g  m'^  increased  seed  number  per  unit  area. 
Fruit  and  seed  numbers  may  be  associated  with  the  availability  of  assimilate 
during  flowering  and  fruit  set  (Stephenson  1981;  Heitholt  et  al.,  1985).    Christy 
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and  Porter  (1982)  further  point  out  that  seed  number  per  unit  area  is  closely 
associated  with  canopy  photosynthesis  during  flowering  and  pod  set. 

At  plant  population  density  D5  yield  per  plant  continues  to  decrease. 
However,  yield  per  unit  area  may  also  begin  to  decrease.  This  decrease  in  yield 
per  area  is  attributed  to  the  traumatic  events  of  lodging  and  barrenness.  Lodging 
and  barrenness  at  high  densities,  are  important  in  explaining  parabolic 
yield-plant  population  density  responses.  At  high  densities,  plants  tend  to  lodge 
because  of  the  failure  of  their  mechanical  axes  as  they  undergo  plastic 
responses  caused  by  increased  interplant  competition.  Apart  from  the  reduced 
yield  that  occurs  when  mechanical  harvesting  is  used  on  lodged  crops,  lodged 
plants  in  comparison  to  upright  plants  have  lower  yields.  This  yield  reduction  has 
been  attributed  to  reduced  light  interception  in  lodged  plants  (Johnston  and 
Pendleton,  1988).  At  very  high  populations  plants  become  barren.  Barrenness 
can  probably  be  attributed  to  the  low  levels  of  light  interception  per  plant.  Pontes 
and  Ohirogge  (1972)  stated  that  these  barren  plants  utilize  light  and  nutrients,  but 
do  not  produce  any  yield.  At  high  densities,  both  yield  per  unit  area  and  yield  per 
plant  decrease  resulting  in  a  parabolic  response.  Leffel  (1961),  however, 
reported  that  natural  lodging  of  soybeans  did  not  significantly  affect  yield  when 
compared  to  unlodged  plants. 

Duncan's  theory  on  vield-plant  population  densitv  responses 

Duncan  (1986)  developed  a  theory  to  explain  yield-plant  population 
density  responses  based  on  two  postulates.  First,  within  soybean  planting 
patterns  there  is  a  range  of  densities  within  which  seed  yield  per  unit  area 
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increases  without  an  increase  in  light  interception  by  the  fully  developed  leaf 
canopy  (this  occurs  at  growth  stage  R5).  Secondly,  within  limits,  the  greater  the 
vegatative  weight  of  a  soybean  plant  during  the  seed  initiation  period  (which 
begins  at  R5)  the  more  seed  it  will  yield,  all  other  conditions  remaining  the  same. 
Duncan  further  argued  that  seed  yield  per  unit  area  was  affected  by  two  dominant 
considerations:  i)  the  fraction  of  total  photosynthetically  active  radiation  (PAR) 
intercepted  by  the  crop  and  ii),  the  efficiency  with  which  the  intercepted  PAR  is 
used  for  seed  production.  Based  on  the  statements  above,  Duncan  defined  three 
Phases  (which  correspond  to  plant  population  density  ranges),  which  are  used  to 
explain  yield-plant  population  density  responses. 

In  Phase  I,  there  is  no  appreciable  competition  between  the  plants  for 
light  ;  light  intercepted  per  plant  and  seed  yield  per  plant  are  constant  and  seed 
yield  per  unit  area  is  influenced  only  by  the  fraction  of  total  light  intercepted. 
Between  Phases  I  and  II,  there  is  an  increasing  degree  of  mutual  shading  as 
plant  population  density  increases.  Light  interception  per  plant  decreases; 
however,  the  total  light  intercepted  (per  area)  increases.  Increases  in  yield  per 
unit  area  in  this  transitional  region  are  due  to  both  increased  light  intercepted  per 
unit  area,  and  increased  efficiency  with  which  intercepted  light  is  used. 

In  Phase  II,  all  of  the  incoming  radiation  is  intercepted  at  full  canopy 
development  (this  occurs  at  growth  stage  R5).  Therefore,  in  Phase  II  light 
intercepted  at  seed  initiation  is  constant,  and  increases  in  seed  yield  per  unit 
area  in  this  phase  are  due  to  the  increased  efficiency  with  which  intercepted  light 
is  used  to  produce  seed.  Phase  II  ends  where  seed  yield  per  unit  area  reaches  a 
maximum. 
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In  Phase  III  seed  yield  per  unit  area  remains  constant  at  a  maximum  for 
the  cultivar  and  the  environment.  Duncan  (1986)  was  able  to  identify  Phases  II 
and  III  using  data  from  Wiggans  (1939)  and  Parks  et  al.  (1983).  Egli  (1988a)  in 
his  work  with  'Cumberland'  found  that  at  R5  and  a  density  of  5.1  plants  m  '^  there 
was  an  apparent  increased  efficiency  with  which  intercepted  light  was  used  for 
seed  production.  This  finding  supports  the  existence  of  Duncan's  (1986)  Phase 
II.  Photosynthetically  active  radiation  intercepted  attained  a  maximum  of  95%, 
after  which  yield  continued  to  increase  with  density.  The  R5  growth  stage  (seed 
initiation)  is  important  in  light  interception  studies  because  it  is  at  this  stage  that 
LAI  attains  a  maximum  (Keller  et  al.,  1970).  Egli's  (1988a)  data  supports 
Duncan's  (1986)  argument  that  plant  population  densities  providing  complete 
insolation  interception  by  the  fully  developed  leaf  canopy,  may  not  be  high 
enough  to  maximize  yield. 


A  basis  for  changes  in  the  efficiencv  with  which  plants  use  Inputs  for  seed 
production 


In  appreciating  the  possible  effect  of  density-mediated  plasticity  on  the 
efficiency  with  which  a  plant  produces  seed,  an  obvious  point  of  departure  is 
Huxley's  law  of  allometry,  which  states  that  the  relative  distribution  of  biomass 
among  the  parts  of  a  plant  depends  on  plant  size  (Willey  and  Heath,  1969).  Kira 
et  al.  (1956)  also  observed  that  the  weight  of  a  plant  part  could  be  related  to  the 
weight  of  the  whole  plant  via 

^    _  i^^h  w  =  whole  plant  weight 

w    =  weight  of  a  plant  part  k  =  a  constant 
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Plant  population  density  has  an  obvious  effect  on  plant  size,  in  addition  to 
other  photomorphogenic  and  plastic  effects.    A  plant's  biomass  is  allocated  to 
vegetative,  photosynthetic,  and  structural  tissues  and  reproductive  organs  and. 
Obviously,  interplant  competition  has  much  potential  for  affecting  the  efficiency 
with  which  inputs  are  used  by  plants  in  seed  production. 

Several  studies  illustrate  that  soybeans  respond  to  reduced  competition 
by  increasing  the  size  of  the  metabolic  sink  sites  (Buttery,  1969;  Johnston  et  al., 
1969;  Lehman  and  Lambert,  1960  )  where  seed  number  and  seed  growth  rate 
are  increased.  At  higher  densities,  however,  the  observed  tendency  for  plants  to 
utilize  intercepted  light  more  efficiently  may  be  of  evolutionary  significance. 
Soybean  plants  may  be  exhibiting  vestiges  of  an  evolutionary  strategy.  Holler 
and  Abrahamson  (1977),  Ogden,  (1974)  and  Thomas  (1974)  reported  that  at  high 
densities  plants  tend  to  maximize  seed  production.  Abrahamson  (1975) 
suggested  that  this  is  a  strategy  to  promote  seed  dispersal  away  from  an 
unfavorable  site.  Egli's  (1988b)  work,  however,  suggested  that  the  partitioning  of 
assimilate  during  this  growth  stage  (flowering  and  fruit  set),  may  not  be  altered  by 
changes  in  plant  size  and  plant  population  density.  Shibles  and  Weber  (1966) 
suggested  that  at  higher  populations,  the  period  of  vegetative  production  is 
longer  than  at  lower  densities.  As  a  consequence,  the  vegetative  production 
period  encroaches  upon  the  seed  production  period,  resulting  in  competition 
within  the  plant  for  available  carbohydrates.  Hence  there  is  less  carbohydrate  for 
seed  production.  Whatever  the  possible  basis  for  changes  in  seed  yield 
efficiency,  it  is  probably  one  of  the  least  understood  areas  in  yield-plant 
population  density  relations. 


CHAPTER  3 
MATERIALS  AND  METHODS 


The  experimental  procedures  used  in  the  field  are  described  prior  to 
discussions  of  data  analysis.  The  discussion  on  data  analysis  involves  the 
development  of  an  analytical  framework  for  the  analysis  of  yield-plant  population 
density  data.  The  materials  and  method  used  with  SOYGRO  will  be  discussed 
along  with  crop  simulation  results  in  chapter  6. 

Field  Experiments 

This  study  was  conducted  at  Gainesville,  Florida  (29°  38')  in  1987  and 
1988.  The  soil  was  classified  as  a  Kendrick  fine  sand  (a  loamy  siliceous 
hyperthermic  family  of  Arenic  Paleudults).  The  previous  crop  consisted  of  small 
grains.  Rainfall  data  were  collected  from  a  nearby  Agronomy  Meteorological 
Station  (Fig.  3-1). 

In  all  the  field  experiments  a  square  planting  arrangement  was  used. 
Wiggans  (1939)  found  that  the  nearer  the  arrangement  of  plants  on  a  given  area 
approaches  a  square  pattern,  the  greater  the  yield  per  unit  area.  Safo-Kantanka 
and  Lawson  (1980)  observed  that  the  number  of  pods  per  plant,  seed  size  and 
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Fig.  3-1 .  Precipitation  profile  over  the  1987and1988  growing 
seasons,  Gainesville,  Fl. 
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yield  per  unit  area  increased  as  a  square  planting  pattern  was  approached.  This 
finding  was  corroborated  by  Parvez  et  al.  (1989)  in  their  work.  Buttery  (1969),  in 
citing  Donald  (1963)  and  Hoiiiday  (1963)  stated  that  it  was  difficult  to 
dennonstrate  any  depressing  effects  of  small  deviations  from  squareness  on  yield. 

Square  planting  is  of  both  analytical  and  production  importance.  A 
constant  square  planting  arrangement  ensures  that  the  differences  observed  in 
plant  population  studies,  resulted  from  plant  population  and  not  plant 
arrangement.  The  square  planting  arrangement  minimizes  inter-plant  competition 
and  maximizes  insolation  interception  per  unit  land  area  (Miura  and  Gemma, 
1986;  Harper,  1983). 

In  1987  the  soybean  cultivars  Kirby  [MG  VIM],  Centennial  [MG  VI]  and 
Forrest  [MG  V]  were  planted  on  27-29  May,  using  planting  boards  which  gave  a 
square  planting  arrangement.  The  inter-row  distances  of  30,  20  and  15  cm  gave 
the  plant  population  densities  of  11.1,  25.0  and  44.4  plants  m"2,  respectively. 
Vermiculite  was  used  to  fill  seeded-hoies  to  ensure  even  germination.  The 
experimental  design  was  a  split-plot  with  plant  population  density  as  the  main 
plot  and  cultivar  as  the  subplot.  Subplots  measured  3.7  by  3.0  m.  The  resulting 
three  main  plots  with  their  combined  total  of  nine  subplots  were  replicated  four 
times.  Fertilizer,  10-10-10  (N-P2O5-K2O)  was  applied  at  a  rate  of  785  kg  ha"""  at 
seed-bed  preparation.  An  overhead  sprinkler  irrigation  system  was  used  during 
periods  of  water  stress.  The  plots  were  hand  weeded.  In  cases  where  gaps 
appeared  in  stands,  soybean  plants  were  transplanted  in  from  border  rows. 
Transplantation  often  resulted  in  weak  plants. 
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The  corn  ear  worm,  Heliothis  zea  (Boddie)  was  controlled  with  lannate, 
(methomyl  S-methyl  N  [(methylcarbamyl)  oxy]  thioacetamidate  which  was  applied 
in  a  solution  containing  2  g  L"''  of  active  ingredient  at  a  rate  of  140  L  ha""".  The 
cyst  nematode,  Heterodera  glycines  affected  crop  productivity,  especially  for  the 
less  nematode  resistant  cultivars,  Forrest  and  Centennial. 

The  following  reproductive  growth  stages  of  soybean  described  by  Fehr  et 
al.  (1971)  were  used  in  sampling  procedures: 

Rl-One  flower  at  any  node 

R2-Flower  at  node  immediately  below  the  uppermost  nodes  with  a 

completely  unrolled  leaf. 
R3-Pod  0.5  cm  long  at  one  of  the  four  uppermost  nodes  with  a 

completely  unrolled  leaf 
R4-Pod  2  cm  long  at  one  of  the  four  uppermost  nodes  with  a 

completely  unrolled  leaf. 
R5-Beans  beginning  to  develop  (can  be  felt  when  the  pod  is  squeeezed 

at  one  of  the  four  uppermost  nodes  with  a  completely  unrolled  leaf. 
R6-Pod  containing  full  size  green  beans  at  one  of  the  four  uppermost 

nodes  with  a  completely  unrolled  leaf. 
R7-Pods  yellowing;  50%  of  leaves  yellow:  physiological  maturity. 
R8-95%  of  pods  brown:  harvest  maturity. 

Whole  plants  were  hand  sampled  at  the  R5,  R6  and  R8  growth  stages  (89, 
132  and  144  d  after  planting,  respectively).  The  samples  were  dried  to  constant 
weight  in  a  forced-draft  oven  at  60  °C.  Total  dry  weights  and  where  appropriate, 
pod  and  seed  weights  and  numbers  were  determined. 

Because  of  the  severe  effects  of  nematode  infestation  on  soybean 
production  in  1987,  the  nematode  resistant  cultivar  Kirby  was  used  in  1988. 
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Plantings  took  place  between  28  -29  June,  with  planting  boards,  which  gave  a 
square  planting  pattern.  The  plant  population  densities  of  2.0,  4.0,  6.3,  11.1, 
16.0,  25.0,  44.4,  and  69.4  plants  m"2  were  achieved.  Seeded-holes  were  filled 
with  builder's  sand  to  ensure  even  germination. 

The  experimental  design  was  a  randomized  complete  block  with  plots 
measuring  4.3  by  3.7  m.  Each  block  consisted  of  eight  plots  and  was  replicated 
four  times.  Fertilizer,  0-10-20  (N-P2O5-K2O)  was  applied  at  a  rate  of  504  kg 
ha"^.  Aldicarb,  [2-methyl-2  (methylthio)  propionaldehyde  o-(methylcarbamyl) 
oxime]  was  applied  at  a  rate  of  27  kg  ha'^ ,  at  seed  bed  preparation  for  nematode 
control.  The  herbicide  alachlor,  2-chloro  2'  6'-diethyl-N-(methoxymethyl) 
acetanilide  was  applied  two  days  after  planting  at  a  rate  of  4.67  L  ha'\  in  a 
solution  containing  480  g  L  '^  of  active  ingredient.  An  overhead  sprinkler 
Irrigation  system  was  used.  Lannate  was  effective  in  controlling  the  corn  ear 
worm,  Heliothis  zea  (Boddie).  however,  it  was  ineffective  in  controlling  the  severe 
infestation  by  white  flies,  Aleyrodidae  aleurocanphus  woglumi.  Ashby.  The  'sooty 
mold',  Capnodium  spp.  also  blackened  leaves  as  they  thrived  on  carbohydrate 
exudates  of  the  white  flies.  It  has  been  shown  that  the  'sooty  mold'  interferes  with 
normal  photosynthesis  (Vaishampayan  and  Kogan,  1980).  Arioglu  et  al.  (1989) 
found  that  there  was  a  significant  negative  correlation  between  white  fly 
infestation  and  seed  yield. 

At  R3  and  R5,  (56  and  68  d  after  planting,  respectively)  light  interception 
above  and  below  the  canopy,  as  well  as  reflectance,  were  measured,  at  or  near 
solar  noon,  with  a  line  quantum  sensor  (LI-COR  model  188).  Whole  plant 
samples  were  taken  at  R3,  R5  and  R8.    It  was  not  possible  to  take  samples 
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between  R5  and  R8  due  to  the  deterioration  of  vegetative  structures  brought  on 
by  white  flies  and  extremely  wet  conditions  (Fig.  3-1). 

Pods  and  seeds  were  retrieved  at  R8.  Subsannples  were  separated  into 
their  component  parts  (stems  and  petioles  were  kept  together).  Leaf  area  was 
determined  using  a  leaf  area  meter  (LI-COR  3100).  Samples  were  dried  in  a 
forced-draft  oven  at  60  °C  to  constant  weight.  Vegetative  dry  weights,  pod  and 
seed  weights  and  numbers  were  measured. 

Subsamples  of  stems  and  leaves  from  R5  were  analyzed  for  N  and  P. 
The  samples  were  chopped  in  a  hammer  mill,  and  then  ground  in  a  Wiley  mill, 
using  a  1  mm  screen.  The  samples  were  digested  using  a  modification  of  the 
aluminum  block  digestion  procedure  of  Gallaher  et  al.  (1975).  To  0.3g  of  sample, 
3.2  g  of  9:1  K2S04:CuS04,  were  added  and  digestion  carried  out  for  4  h  at 
400  °C,  using  10  ml  H2SO4  and  2  ml  H202.  Ammonia  in  the  digestate  and  P 
were  determined  by  semiautomated  colorometery  (Hambleton,  1977)  using  a 
Technicon  Autoanalyzer  II.  The  plant  variables  measured  during  1987  and  1988 
are  provided  in  Table  3-1. 

Data  Analvsis 

The  analytical  framework  developed  here  for  the  analysis  of  yield-plant 
density  data  is  based  on: 

i)  the  identification  of  the  biological  and  ecological  variables 

and  relationships  existing  in  yield-plant  density  studies. 
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Table  3-1.     Summary  of  the  plant  variables  measured  during  1987  and 
1988. 


1 987- — 

yada!2l£& 

Growth  staae 

Total  dry  weight  (g  m"^) 

R5,R6,R8 

Seed  dty  weight  (g  m"^) 

R6,R8 

Pod  dry  weight  (gm^) 

R6,R8 

Seed  number  (m"2) 

R6,R8 

Pod  nurrber  ( m'^) 

R6,R8 

.       .    .           -1QfW-    - 

"        -         -         ■     "        "        — — -        15700 

Total  dry  weight  (g  m"^) 

R3,R5 

Leaf  dry  weight  (g  m"^) 

R3,R5 

Stem  dry  weight  (g  m'^) 

R3,R5 

Leaf  area  (m^) 

R3,R5 

Percent  light  intercepted 

R3,R5 

Leaf  N  concentration  (g  kg"^) 

R5 

Stem  N  concentration  (g  kg'^) 

R5 

Leaf  P  concentration  (g  kg  ''^) 

R5 

Stem  P  concentration  (g  kg  ''^) 

R5 

Seed  dry  weight  (g  m'^) 

R8 

Pod  dry  weight  (g  m'^) 

R5,R8 

Seed  number  (m"^) 

R8 

Pod  number  (m'^) 

R5,R8 
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ij)  the  use  of  analyses  of  variance  (ANOVA),  linear  regression 

and  correlations  to  determine  the  importance  and  strength 
of  these  relationships, 

iii)  the  cross  referencing  of  plant  population  density  effects 

and  yield  associations  to  assess  the  effect  of  plant 
population  density  on  a)  ecological  and  b)  biological 
relationships  which  determine  or  are  associated  with  yield. 

The  basic  data  set   (Table  3-1.)  was  expanded  to  include  the  following 
calculated  variables: 

Stem  fraction  =  Stem  weight  /  Total  dry  weight 

Leaf  fraction  =  Leaf  dry  weight  /  Total  dry  weight 

Leaf  area  index(LAI)  =  Leaf  area  /  Harvested  area 

Specific  leaf  area  (m^  g-l )  =  Leaf  area  /  Leaf  dry  weight 

Specific  leaf  weight  (g  m"^)  =  Leaf  dry  weight  /  Leaf  area 

Seed  number  Pod'^  =  Seed  number  /  Pod  number 

Weight  seed'""  (g)  =  Seed  dry  weight  /  Seed  number 

Leaf  N  (g  m'^)  =  Leaf  dry  weight  X  Leaf  N  concentration  X  0.01 

Leaf  N  (g  m-2(leaf  area))  =  Leaf  N  (g  m  "2)  /  LAI 

Stem  N  (g  m"2)  =  Stem  dry  weight  X  Stem  N  concentration  X  0.01 
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Leaf  P  (g  m'^)  =  Leaf  dry  weight  X  Leaf  P  concentration  X  0.01 

Leaf  P  (g  m-2(ieaf  area))         =  Leaf  P  (g  m'^)  /  LAI 

Stem  P  (g  m"2)  =  Stem  dry  weight  X  Stem  P  concentration  X  0.01 


The  Evaluation  of  Variables  and  Relationships  Operating  in  Yield-Plant 
Population  Density  Studies 


In  studying  the  yield  determining  relationships  (both  ecological  and 
biological)  existing  in  crops  which  are  subject  to  varying  plant  population 
densities,  these  relationships  first  need  to  be  classified;  secondly,  they  need  to 
be  measured;  and  thirdly,  their  influence  on  yield  needs  to  be  evaluated. 

The  classification  of  variables  and  relationships 

In  investigating  the  process  of  yield  achievement,  variables  and  the 
relationships  between  variables,  may  contribute  to  or  affect  three  spheres  of  crop 
productivity.  These  spheres  are,  photosynthetic  activity,  nutrient  remobilization 
and  reproductive  parameters  (i.e.,  the  size  and  number  of  'sinks').  Variables  and 
their  inter-relationships  in  contributing  to  these  various  productive  spheres,  in 
yield-plant  density  studies,  may  be  either  plastic  or  non-plastic.  For  example, 
weight  seed'^  is  known  to  be  non-plastic,  being  largely  under  genetic  and  not 
under  environmental  control.  In  this  investigation,  variables  and  relationships 
described  as  not  being  related  nor  associated  with  yield,  refer  specifically  to  the 
tendency  of  these  variables  not  to  undergo  changes  similar  in  magnitude  to 
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changes  in  the  magnitude  of  yield.   The  extent  of  the  contribution  of  'plastic'  and 
'non-plastic'  variables  to  yield  is  not  being  investigated. 

The  number  of  relationships  between  measured  responses  (including 
those  which  have  been  calculated)  is  given  by, 

nCr    =      nPrC!)"''       =       n!  (r!(n-r)!)-1 

pC^   =      the  total  number  of  relationships,  taking  two  responses  at  a 

time 
pPp    =      the  total  number  of  permutations  of  n  responses  taken  r  at  a 

time 
n       =      number  of  measured  or  calculated  responses 
r        =      the  number  of  responses  assessed  at  a  time  (i.e.,  two) 

For  the  1988  experiment,  there  is  a  combined  total  of  28  variables  from 
the  R5  and  R8  stages.  The  total  number  of  possible  relationships  between  the  R5 
and  R8  variables  is  378  (  28!  (2!  x  26!)  ■■•).  There  are  406  (378  +  28  (plant 
population  density  relationships))  relationships  which  may  be  analyzed. 

Techniques  for  determining  relationships  in  vield-plant  densitv  studies 

The  three  basic  techniques  which  will  be  used  in  measuring  the  extent  of 
the  relationships  operating  in  yield-plant  density  studies  are  the  analysis  of 
variance  (ANOVA),  linear  regression  and  correlations.  The  ANOVA  will  be  used 
to  determine  density  and  cultivar  effects  on  measured  variables.  Differences 
between  cultivars  will  be  determined  using  Duncan's  Multiple  Range  Test 
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(DMRT).  Ordinarily,  It  is  the  aim  of  linear  regression  analysis  to  describe  a 
functional  relationship  between  Y  and  X  (i.e.,  plant  population  density)  (  Sachs, 
1982).  However,  results  of  the  linear  regression  analyses  obtained  in  this 
investigation  will  not  be  seen  as  describing  the  functional  relationship  between 
plant  population  density  (X)  and  a  given  variable  (Y).  This  is  because  linear 
regressions  are  not  always  appropriate  for  providing  exact  relationships 
between  variables  in  complex  biological  systems.  Linear  regression  analyses  of 
plant  population  density  effects,  are  used  here  strictly  as  statistical  tests,  which 
establish  the  extent  of  the  impact  of  plant  population  density  on  a  given  variable. 
Therefore,  of  specific  interest  are  the  R^  values  obtained  in  linear  regressions. 
In  using  linear  regression  as  a  statistical  test,  logarithmic  transformations  serve 
the  dual  function  of  linearizing  the  data  and  stabilizing  the  variance  (Montgomery, 
1984).  The  Pearson  product  moment  correlation  coefficient  (r)  which  will  be 
used  here  as  a  statistic  for  the  strength  of  the  relationship  between  two  variables. 
Correlations  will  be  used  to  assess  the  strength  and  nature  (whether  it  is  positive 
or  negative)  of  primary,  secondary  and  tertiary  relationships. 

The  evaluation  of  the  relationships  operating  in  vield-plant  density  studies 

The  evaluation  of  the  ecological  and  biological  relationships  which  may 
be  yield-determining  will  be  conducted  by  cross  referencing  i),  plant  population 
density  effects  on  variables,  and  ii),  the  association  of  these  variables  with  yield. 
The  first  step  is  the  determination  of  the  strength  of  the  effect  of  plant  population 
density  on  measured  variables.     Secondly,  the  strength  of  the  association 
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between  measured  variables  and  yield  per  unit  area  is  determined.  These 
relationships  are  then  categorized  based  on  R^  (from  linear  regression)  and  r 
(from  correlations)  values  into  high  (H)  (0.75  -  1.00),  moderate  (M)  (0.50  -  0.74) 
and  low  (L)  (0.00  -  0.49)  categories.  The  next  step  which  will  be  illustrated  using 
only  the  "high"  (H)  and  'low'  (L)  categories,  involves  the  correlation  of  'responses 
classified  with  respect  to  their  relationship  with  yield'  against  'responses 
classified  with  respect  to  their  relationship  with  plant  population  density".  The 
resulting  r,  is  then  used  to  determine  the  possible  contribution  of  primary, 
secondary  and  tertiary  relationships  to  yield  per  unit  area,  in  plant  population 
density  studies. 

Based  on  R^  and  r  values,  variables  can  be  classified  into  those;  affected 
by  plant  population  density  (HDV);  not  affected  by  plant  population  density  (LDV); 
associated  with  yield  (HYV);  and  not  associated  with  yield  (LYV).  In  some 
instances,  HDVs  and  HYVs  may  be  the  same.  The  types  of  relationships  which 
may  exist  between  these  variables  are, 


HDVs  X  HYVs 
HDVs  X  LYVs 
LDVs  X  HYVs 
LDVs  X  LYVs 

The  interpretation  of  the  strength  of  these  relationships  as  determined  by  r 
values  are  outlined  in  Table  3-2.  The  relationships  which  are  of  greatest  interest 
are  those  involving  HDVs  X  HYVs,  and  LDVs  X  LYVs  with  high  r  values.  The 
former  identifies  variables  and  relationships  which  are  particularly  plastic,  and 
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Table  3-2.  A  table  for  the  interpretation  of  plant  population  density  -  yield, 
cross-referenced  data. 


Type  of  relation 


A  high  r 


A  low  r 


HDVs  X  HYVs 


A  strong  relationship  which 
is  most  likely  to  be  yield 
determining  and  which  is 
under  plant  population 
density  effect. 


A  weak  relationship,  however 
the  HYVs  are  greatly  associated 
with  yield,  and  the  HDVs  are 
under  plant  population  density 
effect.  This  relationship  cannot 
be  yield  determining  in  plant 
population  density  studies. 


HDVs  X  LYVs 


A  strong  relationship, 
however,  this  relationship, 
cannot  explain  yield  as  a 
plant  density-mediated 
process. 


A  weak  relationship  not 
capable  of  explaining  yield, 
however,  the  HDVs  are 
under  plant  population  density 
effect. 


LDVs  X  HYVs 


A  strong  relationship, 
however,  this  relationship, 
cannot  explain  yield  as  a 
plant  density-mediated 
process. 


Although  the  HYVs  are 
closely  associated  with  yield, 
this  weak  relationship  is  not 
capable  of  explaining  yield 
as  a  plant  density-mediated 
process. 


LDVs  X  LYVs 


A  strong  relationship  between 
variables  which  are  not 
associated  with  yield  and 
are  not  affected  by  plant 
density. 


A  weak  relationship  between 
variables  which  are  not 
associated  with  yield 
and  are  not  affected 
by  plant  population  density. 
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which  are  most  likely  to  be  responsible  for  the  changes  in  the  magnitude  of  yield 
per  unit  area  with  changes  with  plant  population  density.  The  latter  identifies 
variables  and  relationships,  which  are  not  particularly  plastic  (which  are  probably 
under  much  genetic  control)  and  cannot  explain  changes  in  the  magnitude  of 
yield  per  unit  area  with  changes  in  plant  population  density. 


CHAPTER  4 
RESULTS  AND  DISCUSSION  -1987  EXPERIMENT 


Plant  population  had  no  significant  effect  on  any  of  the  variables 
measured.  In  the  absence  of  significant  plant  population  effect  on  crop  variables, 
the  cross  referencing  of  the  effect  of  plant  population  density  on  variables  with  the 
'yield  association'  of  variables  is  not  possible.  However,  the  cultivar  effect  on  crop 
variables  and  the  effect  of  plant  population  density  at  levels  above  the  5% 
significance  level,  allows  for  the  use  of  ANOVA  and  correlation  analyses  for  yield 
investigation. 

ANOVA  Results 

Plant  Population  Effect 

Plant  population  had  no  effect  on  any  of  the  variables  investigated  at  R5 
(Table  4-1),  R6  (Table  4-3),  and  R8  (Table  4-5).  In  instances  where  this  result 
was  unexpected  (such  as  with  seed  and  total  dry  weight  per  unit  area),  it  may  be 
attributed  to  a  combination  of  i),  the  plant  populations  used  were  not  sufficiently 
discriminating,  and  ii),  nematode  infestation.  The  nematode  resistant  soybean 
varieties  recommended  for  Florida  (IFAS,  1983),  does  not  include  Kirby  and 
Forrest.  Centennial,  however,  is  recommended.  Analyses  conducted  on  root 
and  soil  samples  showed  that  the  cyst  nematode  (Heterodera)  (race  unidentified, 
and   mostly  juveniles),  and   root  knot  nematode   (Meloidogvne  incognita. 
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Table  4-1 .  The  effect  of  plant  population  and  cultivar  on  total  dry 
weight  (g  m"2)  at  R5, 1987. 


Source  of  variation 


Population  5.56 

Cultivar  31.8" 

Population  X  Cultivar  1.96 


*,  **,  significance  at  the  5  and  1%  levels,  respectively. 


Table  4-2  Cultivar  effect  on  total  dry  weight  at  R5, 
1987. 


Cultivar  Total  dry  weight  (g  m'^) 


Kirby  957a* 

Centennial  723b 

Forrest  551c 


*  observations  followed  by  the  same  letter  are  not 
significantly  different  (DMRT,  5%). 
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Table  4-3  The  effect  of  plant  population  and  cultivar  on  soybean  plant  traits  at  R6 
stage,  1987. 


Variable 

Source  of  variation 

F 

Total  dry  weight  (g  m"^) 

Population 

Cultivar 

Population  X  cultivar 

0.21 
10.84" 
1.91 

Pod  weight  (g  m"^) 

Population 

Cultivar 

Population  X  cultivar 

0.15 

12.06" 

2.28 

Seed  dry  weight  (g  m"^) 

Population 

Cultivar 

Population  X  cultivar 

0.27 

10.26" 
1.92 

Pod  number  (m"^) 

Population 

Cultivar 

Population  X  cultivar 

0.34 
1.90 
2.18 

Seed  number  (m"^) 

Population 

Cultivar 

Population  X  cultivar 

0.28 
2.69 
2.24 

Seed  number  Pod "'' 

Population 

Cultivar 

Population  X  cultivar 

2.32 
0.56 
1.27 

Weight  seed  ■'' (g) 

Population 

Cultivar 

Population  X  cultivar 

1.07 
33.37" 
0.73 

*,  ",  significance  at  the  5  and  1%  levels,  respectively. 
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Table  4-4.  Cultivar  effect  on  soybean  plant  traits  at  R6  stage,  1987. 


Cultivar 

Variable 

Kirby 

Centennial 

Forrest 

Total  dry  weight  (g  m'^) 

644  a* 

598  a 

443  b 

Pod  weight  (g  m"^) 

298  a 

318  a 

223  b 

Seed  dry  weight  (g  m"^) 

159  b 

198  a 

141  b 

Pod  number  (m"^) 

1236  a 

1213a 

1056  a 

Seed  number  (m"^) 

2400  a 

2294  a 

2015  a 

Seed  number  Pod  '^ 

1.95  a 

1.90  a 

1.94  a 

Weight  seed  ""'(g) 

0.07  b 

0.09  a 

0.07  b 

*  observations  followed  by  the  same  letter  are  not  significantly 
different  (DMRT,  5%). 
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Table  4-5.  The  effect  of  plant  population  and  cultivar  on  soybean  plant  traits  at  R8 
stage,  1987. 


Variable 

Source  of  variation 

F 

Pod  weight  (g  m'^) 

Population 
Cultivar 

2.14 
11.82" 

Population  X  cultivar 

1.92 

Seed  dry  weight  (g  m"^) 

Population 
Cultivar 

1.08 
12.76" 

Population  X  cultivar 

1.85 

Pod  number  (m*^) 

Population 
Cultivar 

0.73 
1.95 

Population  X  cultivar 

4.03* 

Seed  number  (m'^) 

Population 
Cultivar 

1.53 
2.29 

Population  X  cultivar 

4.36* 

Seed  number  Pod  "^ 

Population 
Cultivar 

0.57 
8.67** 

Population  X  cultivar 

0.42 

Weight  seed  ""I  (g) 

Population 
Cultivar 

0.40 
14.69" 

Population  X  cultivar 

1.73 

',  **,  significance  at  the  5  and  1%  levels,  respectively. 
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'Southern')    populations  were  sufficient  to  hinder  grain  production  in  all  three 
cultivars. 

At  R5,  total  dry  weight  (g  m  '^)  tended  to  increase  with  population  (Fig. 
4-1).  Although  such  an  observation  is  consistent  with  the  results  of  Udoguchi  and 
McCloud  (1986),  Blumenthal  et  al.  (1988)  and  Parvez  et  al.  (1989),  the 
population  effect  was  not  significant  (p  =  0.07)  (Table  4-1).  There  was  no 
significant  difference  in  seed  dry  weight  between  plant  populations  at  R8  (Fig. 
4-2)  (Table  4-5),  and  this  can  be  attributed  to  the  effect  of  plant  population  density 
on  yield  components.  The  number  of  seeds  per  pod  (R8)  was  not  affected  by 
plant  population  (Fig.  4-3)  (Table  4-5).  This  result  corroborates  those  of 
Chaudhry  and  Cheema  (1985)  and  Dominguez  and  Hume  (1978).  Weight  per 
seed  (R8)  (g)  remained  unchanged  with  plant  population  (Fig.  4-4)  (Table  4-5), 
and  this  is  similar  to  the  results  obtained  by  Dominguez  and  Hume  (1978),  and 
Fontes  and  Ohirogge  (1972).  Pod  number  m"^  was  unaffected  by  plant 
population  density  (Fig  4-5)  (Table  4-5).  With  none  of  the  yield  components 
being  affected  by  plant  population,  the  absence  of  plant  population  effect  on  seed 
dry  weight  (g  m"^)  at  R8  is  expected. 

Cultivar  Effect 

Total  dry  weight  (g  m"^)  at  R5  was  affected  by  cultivar  (Table  4-1)  with 
Kirby,  Centennial  and  Forrest  providing  970,  723  and  551  g  m'^,  respectively 
(Table  4-2  ).  At  R6,  pod  number  m'^,  seed  number  m"^  and  seed  number  pod"^ 
were  not  affected  by  cultivar  (Table  4-3).  However,  total  dry  weight  (g  m'^),  pod 
dry  weight  (g  m"^),  seed  dry  weight  (g  m'^)  and  weight  per  seed  (g),  were  all 
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Fig.  4-1.  The  effect  of  plant  population  density  on  total  dry 
weight  at  the  R5  stage  in  soybean,  1 987. 
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Fig.  4-2.  The  effect  of  plant  population  density  on  seed  dry 
weight  of  soybean  at  maturity,  1 987. 
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soybean  at  maturity,  1 987. 
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Fig.  4-4.  Theeffectofplantpopulationdensity  on  weight  seed     of 
soybean  at  maturity,  1 987. 
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affected  by  cultivar.  These  results  are  summarized  in  Tables  4-3  and  4-4.  At  R8 
(Table  4-5),  all  variables  were  affected  by  cultivar.  The  interaction  of  cultivar  and 
population  on  seed  number  m'^,  and  pod  number  m"^  warranted  further 
investigation;  upon  analysing  the  cultivars  individually,  for  both  pod  number  m"^ 
(Table  4-6)  and  seed  number  m"^  (Table  4-7)  there  was  no  population  effect  in 
either  case.  At  each  plant  population,  there  was  no  cultivar  effect  on  pod  number 
m"2  (Table  4-8)  and  seed  number  m'^  (Table  4-9).  A  summary  of  the  cultivar 
effects  at  RB  is  provided  in  Table  4-10. 

An  Explanation  for  Cultivar  Yield  Differences  Based  on  ANOVA 

In  Table  4-10.  a  comparison  of  seed  dry  weight  (RB)  (g  m"^)  between 
Kirby,  Centennial  and  Forrest  shows  that  yields  were  26B,  270  ,  and  17B  g  m  '2, 
respectively,  with  yield  of  Forrest  being  significantly  lower.  A  review  of  yield 
components  and  total  dry  weight  (R5)  (g  m'2)  shows  that  pod  number  m"^,  and 
seed  number  pod  '\  are  unable  to  account  for  differences  in  seed  dry  weight  (g 
m"2).  Weight  seed' ""(g)  and  total  dry  weight  (R5)  (g  m"^)  appear  to  be 
responsible  for  the  seed  dry  weight  (g  m"^)  differences  observed  for  these 
genotypes. 

It  has  been  shown  that  pod  number  m'^,  and  seed  number  pod"^,  are 
unable  to  account  for  differences  in  seed  dry  weight  (g  m'^).  Kirby  and 
Centennial  are  from  higher  maturity  groups  than  Forrest  and,  therefore,  have 
longer  grain  filling  periods,  and  this  may  have  attributed  to  Kirby  and  Centennial 
having  larger  seeds.  Also,  of  the  three  cultivars,  Forrest  was  the  least  nematode 
resistant  and  this  may  have  taken  away  further  from  the  ability  of  Forrest  to 


53 


Table  4-6.  The  effect  of  plant  population  (X)  on  pod  number  (m'^)  (Y)  for  each 
cultivar  at  R8  stage,  1987. 


Cultlvar  Statistical  model  R^ 


Kirby  InY  =  0.01  InX  +  7.29  0.00 

Centennial  InY  =  -0.32  InX  +  8.27  0.35 

Forrest  InY  =  0.26  InX  +  6.29  0.37 


',**  significance  at  the  5  and  1%  levels,  respectively. 


Table  4-7.  The  effect  of  plant  population  on  seed  number  (m"^)  for  each  cultivar 
at  R8  stage,  1987. 

Cultivar  Statistical  model  R^ 


Kirby  In  Y  =  -0.05  InX  +  8.03  0.02 

Centennial  In  Y  =  -0.32  InX  +  8.94  0.37 

Forrest  In  Y  =  0.23  InX  +  7.04  0.35 


',**  significance  at  the  5  and  1%  levels,  respectively. 
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Table  4-8.  The  effect  of  cultivar  on  pod  number  (m"^)  for  each  plant  population  at 
R8  stage, 1987. 


Plant  population  Source  of  variation  F 

44.4  plants  m'2  Cultivar  3.53 

25.0  plants  m"2  Cultivar  2.20 

11.1  plants  m"2  Cultivar  4.22 
',**  significance  at  the  5  and  1%  levels,  respectively. 


Table  4-9.  The  effect  of  cultivar  on  seed  number  (m'^)  for  each  plant  population 
at  R8  stage,  1987. 


Plant  population           Source  of  variation  F 

44.4  plants  m"2            Cultivar  3.01 

25.0  plants  m"2            Cultivar  3.11 

11.1  plants  m"2           Cultivar  4.50 

*,**  significance  at  the  5  and  1%  levels,  respectively. 
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Table  4-10.  Cultivar  effect  on  soybean  plant  traits  at  R8  stage,  1987. 


Variable 


Cultivar 


Kirby 

Centennial 

Forrest 

405  a 

385  a 

266  b 

268  a 

270  a 

178  b 

1493  a 

1499  a 

1265  a 

2701  a 

2964  a 

2457  a 

1.81  b 

1.99  a 

1.95  a 

0.10  a 

0.09  a 

0.07  b 

Pod  weight  (g  m'^) 
Seed  dry  weight  (g  m"^) 
Pod  number  (m"^) 
Seed  number  (m"2) 
Seed  number  Pod  "^ 
Weight  seed"^(g) 


*  observations  followed  by  the  same  letter  are  not  significantly  different 
(DMRT.  5%) 


56 


perform.  Forrest  also  had  the  lowest  vegetative  dry  weight  at  R5,  and  this  means 
that  in  comparison  to  Kirby  and  Centennial,  less  nutrients  were  available  for 
remobilization. 

Correlation  Analvsis 

The  correlation  of  variables  at  R5  and  R6,  and  R5  and  R8,  with  seed  dry 
weight  (g  m"^)  and  yield  components,  are  provided  in  Tables  4-11  and  4-12, 
respectively.  Total  dry  weight  (g  m"^)  at  R5  was  poorly  correlated  with  seed  dry 
weight  (g  m'^)  at  R6  (r  =  0.34,  p  <  .05).  Total  dry  weight  (g  m"^)  at  R5  was 
moderately  correlated  with  seed  dry  weight  (g  m"^)  at  R8  (r  =  0.56,  p  <  .01 ).  Total 
dry  weight  (g  m'^)  at  R5  was  also  moderately  correlated  with  all  yield 
components  at  R8. 

At  both  R6  and  R8,  pod  dry  weight  (g  m"^),  pod  number  (m'2)  and  seed 
number  (m'^)  were  highly  correlated  with  seed  dry  weight  (g  m'^). 
Basuchaudhari  (1987),  Kaw  and  Menon  (1972),  Pontes  and  Ohirogge  (1972), 
and  Parvez  et  al.  (1989)  had  similar  findings.  Weight  seed  "^  (g)  was 
moderately  correlated  with  seed  dry  weight  (g  m'^)  at  both  R6  and  R8.  These 
results  are  contrary  to  those  obtained  by  Board  (1985),  Tandojam  (1986), 
Johnson  et  al.  (1955).  However,  Moorthy  (1952)  found  that  there  was  a  high 
correlation  between  weight  seed  "^(g)  and  seed  dry  weight  (g  m"^).  At  R6  there 
was  moderate  correlation  between  seed  dry  weight  (g  m"^)  and  seed  number 
pod*^,  this  result  was  also  obtained  by  Ma  (1946).  At  R8,  seed  number  pod"^ 
was  not  correlated  with  seed  dry  weight  (g  m"^).  This  result  corroborates  those 
obtained  by  Hartwig  and  Edwards  (1970),  Board  (1985),  Tandojam  (1986)  and 
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Table  4-1 1 .  Correlation  coefficients  for  seed  yield  and  its  fundamental 
components  (XI)  at  R6  with  variables  (X2)  at  R5  and  R6,  1987. 

XI  X2  r 

Seed  dry  weight  (R6)(g  m"^)  Total  dry  weight  (R5)(g  m"^)  0.34* 

Total  dry  weight  (R6)(g  m"^)  0.82" 

Pod  dry  weight  (R6)(g  m"^)  0.93** 

Pod  number  (R6)(m"^)  0.83** 

Seed  number  (R6)  (m"^)  0.81  ** 

Seed  number  Pod"^  (R6)  -0.48* 

Weight  seed  '■"  (R6)  (g)  -0.56** 

Pod  number  (R6)  (m'^)  Total  dry  weight  (R5)(g  m"^)  0.55** 

Total  dry  weight  (R6)(g  m"^)  0.87** 

Pod  dry  weight  (R6)  (g  m'^i  0.88** 

Seed  dry  weight  (R6)  (g  m"^)  0.83** 

Seed  number  (R6)  (m"'^)  0.97** 

Seed  number  Pod" '  (R6)  -0.56** 

Weight  seed  "'' (R6)(g)  0.06 

Seed  number  Pod"""  (R6)  Total  dry  weight  (R5)(g  m"^)  -0.21 

Total  dry  weight  (R6)(g  m"^)  -0.35 

Pod  dry  weight  (R6)  (g  m'^i  -0.42* 

Seed  dry  weight  (R6)  (g  m'^)  -0.48* 

Pod  number  (R6)  (m'^)  -0.56** 

Seed  number  (R6)  (m'^)  -0.34 

Weight  seed  "''  (R6)  (g)  -0.35 

Weight  seed  "''  (R6)  (g)  Total  dry  weight  (R5)(g  m"^)  -0.12 

Total  dry  weight  (R6)(g  m"^)  0.14 

Pod  dry  weight  (R6)  (g  m'^j  0.33 

Seed  dry  weight  (R6)  (g  m'^)  0.56** 

Pod  number  (R6)  (m"^)  0.06 

Seed  number  (R6)  (m'^)  -0.03 

Seed  number  (R6)  Pod"''  -0.35 

*,**  significance  at  the  5  and  1%  levels,  respectively. 
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Table  4-1 2.  Correlation  coefficients  for  the  fundamental  yield  components  (X1 )  at 
stage  R8  with  variables  (X2)  at  R5  and  R8, 1 987. 


X1 


X2 


Seed  dry  weight  (R8)(g  m"*^) 


Pod  number  (R8)  (m"^) 


Seed  number  Pod'"'  (R8) 


Weight  seed  '"•  (RB)  (g) 


Total  dry  weight  (R5)  (g  m'^) 
Pod  dry  weight  (R8)  (g  m'^) 
Pod  number  (R8)  (n\'\ 
Seed  number  (R8)  (m'^) 
Seed  number  Pod"^  (R8) 
Weight  seed "''  (R8)  (g) 

Total  dry  weight  (R5)  (g  m"^) 
Pod  dry  weight  (R8)  (g  m'^^ 
Seed  dry  weight  (R8)  (g  m'^) 
Seed  number  (R8)  (m''') 
Seed  number  Pod'^  (R8) 
Weight  seed '''  (R8)  (g) 

Total  dry  weight  (R5)  (g  m"^) 
Pod  dry  weight  (R8)  (g  m'\ 
Seed  dry  weight  (R8)  (g  m'^) 
Pod  number  (R8)  (m'^) 
Seed  number  ms)  (m"^) 
Weight  seed  "''  (R8)  (g) 

Total  dry  weight  (R5)  (g  m"^) 
Pod  dry  weight  (R8)  (g  m"^^ 
Seed  dry  weight  (R8)  (g  m'^) 
Pod  number  (R8)  (m'\ 
Seed  number  (R8)  (m"^) 
Seed  number  Pod" '  (R8) 


0.56" 
0.98" 
0.83" 
0.83" 
-0.13 
0.45* 

0.43* 
0.79** 
0.83** 
0.97** 
-0.28 
-0.06 

-0.56* 

-0.15 

-0.13 

-0.28 

-0.04 

-0.16 

0.55** 

0.47* 

0.45* 

-0.06 

-0.11 

-0.16 


*,**  significance  at  the  5  and  1%  levels,  respectively. 
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Johnson  et  al.  (1955).  Among  yield  components  at  both  R6  and  R8,  seed  number 
pod"^  and  weight  seed"^(g)  are  the  least  correlated  with  the  variables  measured, 
and  this  may  indicate  that  these  yield  components  are  largely  under  genetic 
control. 


An  Explanation  for  Cultivar  Yield  Differences  Based  on  Correlation 
Analyses 


Based  on  DMRT,  it  has  already  been  shown  that  weight  seed  "^(g)  and 
total  dry  weight  (R5)  (g  m"^)  appear  to  be  associated  with  cultivar  seed  dry  weight 
(g  m'2)  differences  (Table  4-11).  However,  a  number  of  important  points  are 
brought  out  in  attempting  to  explain  differences  observed  in  seed  dry  weight 
(g  m'2)  between  Kirby,  Centennial  and  Forrest  using  correlations.  First,  weight 
seed  "''  and  total  dry  weight  (R5)  (g  m"2),  which,  based  on  DMRT,  can  explain 
cultivar  differences  in  seed  dry  weight  (g  m"^)  are  not  highly  correlated  with  seed 
dry  weight  (g  m'2)  (Table  4-13).  This  raises  the  point  that  the  investigation  of  the 
effect  of  individual  variables  on  yield,  as  i),  with  DMRT  analyses  (Table  4-10)  or 
ii),  the  correlation  of  a  single  variable  with  seed  dry  weight  (g  m"^)  has  limitations. 
Cultivar  differences  observed  in  seed  dry  weight  (g  m"^)  are  evidently  attributable 
to  a  number  of  variables  acting  together,  and  the  investigation  of  the  effect  of  one 
variable  on  seed  dry  weight  (g  m"^)  yields  an  incomplete  result. 

Secondly,  in  Table  4-13,  the  combined  data  from  all  cultivars,  when 
regressed  against  seed  dry  weight  (g  m'^),  provides  r  values  which  do  not 
represent  r  values  obtained  for  each  cultivar.  This  discrepancy  is  particularly  the 
case  for  seed  pod  "^  and  total  dry  weight  (g  m"^).  Therefore,  the  individual  yield 
dynamics  of  cultivars  and  the  differential  response  of  cultivars  to  plant  population 
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Table  4-1 3.  Mean  values  and  correlation  coefficients  for  R5  total  dry  weight 
(g  m*2),  seed  dry  weight  (g  m'^)  and  its  fundamental  components  at  R8,1987. 


Cuttvar 


Variable 


jsm. 


Centennial 


Forrest 


All 
cultivars 


Seed  dry  weighty  m"2)  268a+  270a 

Pod  number  (m-2)  1493a  (.77*)@  1499a  (.96**) 

Seed  Pod  ■■•  1 .81  b  (.33)  1 .99a  (-.52) 

Weight  seed -■'(g)  0.10a  (.48)  0.09a  (-.47) 

Total  dry  weight(R5)(g  m-2)  957a  (.30)  723b  (.35) 


178b 

1265a  (.93**) 
1.95a  (-.33) 
0.07b  (-.21) 
551c  (.56) 


(.83**) 
(-.13) 
(.45*) 
(.56**) 


■''  observations  followed  by  the  same  letter  are  not  significantly  different 
(DMRT,  5%). 

^Correlation  coefficients  (r)  of  variables  with  seed  dry  weight  (g  m'^)  are  in 
parentheses. 

*,**  significance  at  the  5  and  1%  levels,  respectively. 


Table  4-14.  Statistical  models  for  per  plant  responses  (Y)  of  soybean  with  plant 
density  (X) 


Variable  (Y) 


Statistical  model 


r2 


Vegetative  dry  weight  Plant  "^  (R5) 
Pod  dry  weight  Plant  ''^  (R8) 
Seed  dry  weight  Plant  •"■  (RB) 
Pod  number  Plant  ""I  (RB) 
Seed  number  Plant  """(RB) 


lnY  = -0.81  lnX  + 5.99 


lnY  =  -1.02lnX  +  5.90 


lnY  =  -1.04lnX-h5.55 


lnY  =  -1.02lnX  +  7.2B 


lnY  =  -1.04lnX  +  8.01 


0.74" 
0.84" 
0.80" 
0.82** 
0.84" 


significance  at  the  5  and  1%  levels,  respectively. 
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density,  indicates  that  in  depth  yield-plant  population  density  studies  should 
focus  on  a  single  cultivar. 

The  Yield-Plant  Population  Density  Response. 

Seed  weight  (g  m"^)  was  not  affected  by  plant  population  density,  and  no 
lodging  was  observed  at  the  highest  plant  population  densities.  However,  the 
vegetative  dry  weight  (R5)  (Fig  4-6),  seed  weight  and  pod  weight  (Fig.  4-7  ),  and 
seed  and  pod  number  (Fig.  4-8),  on  a  per  plant  basis  were  all  significantly 
affected  by  plant  population  density  (Table  4-14).  The  reduction  in  vegetative  dry 
weight  per  plant  with  density  (Fig.  4-6),  indicates  a  reduction  in  the  amount  of 
nutrients  which  may  be  remobilized  for  seed  production.  Individual  plants  clearly 
adjusted  plastically  to  increased  plant  competition;  at  higher  plant  population 
densities  they  reduced  the  number  of  sink  sites  which  is  consistent  with  reduced 
inputs  at  higher  plant  population  densities.  It  apppeared  that  interplant 
competition  had  Its  greatest  effect  between  the  plant  population  densities  of  11.1 
and  25.0  plants  m"^,  where  vegetative  dry  weight  and  yield  reduction  occurred  at 
the  greatest  rate.  As  population  density  approached  44.4  plants  m'^,  plastic 
adjustments  within  the  plant  reduced  the  rate  of  yield  decrease  with  plant 
population  density-that  is,  the  plants  were  adapting  and  responding  positively  to 
the  negative  forces  of  interplant  competition. 
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Fig.  4-6.  The  effect  of  plant  population  density  on  dry  weight  per 
plant  of  soybean  at  the  R5  stage,  1 987. 
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Fig.  4-7.  The  effect  of  plant  population  density  on  seed  and  pod 
weight  per  plant  of  soybean  at  maturity,  1 987. 
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Fig.  4-8.  The  effect  of  plant  population  density  on  seed  and  pod 
number  per  plant  of  soybean  at  maturity,  1 987. 
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Conclusion 

Although  plant  population  had  no  significant  effect  on  the  variables 
measured  at  R5,  R6  and  R8,  the  effect  of  cultivar  on  variables  allowed  for  yield 
investigations.  Total  dry  weight  (g  m  '^)  was  moderately  correlated  with  all  yield 
components,  however,  it  was  negatively  correlated  with  seed  number  pod  '^. 
Correlation  analyses  on  the  data  sets  for  each  cultivar,  and  for  the  combined 
data  set  of  all  cultivars,  at  both  R6  and  R8,  showed  that  pod  dry  weight 
( g  m  "2),  seed  number  (m"^),  and  pod  number  (m  "2)  were  all  highly  correlated 
with  seed  dry  weight  (g  m'^).  Seed  number  po6'\  was  unaffected  by  plant 
population  density,  and  Kirby  had  the  lowest  value.  However  the  greater  weight 
seed  "''  (g)  of  Kirby,  contributed  to  Kirby  having  a  greater  seed  dry  weight 
(g  m  '2)  at  R8  than  Forrest. 

Weight  seed  '^  (g)  was  neither  plastic  nor  was  it  highly  correlated  with 
seed  dry  weight  (g  m'^).  However,  among  yield  components,  weight  seed  ' ' 
accounts  for  the  yield  difference  observed  in  cultivar.  This  apparent  contradiction 
in  the  influence  of  weight  seed  '^  (g)  highlights  the  following  point.  Weight 
seed'^  (g)  determines  final  seed  dry  weight  (g  m"^);  however,  being  under 
genetic  control,  it  can  account  for  differences  in  genotypes,  but  weight  seed  '^ 
(g)  is  unable  to  account  for  the  plastic  response  of  seed  weight  (g  m"^). 

The  combination  and  subsequent  analysis  of  data  obtained  from  different 
cultivars  in  yield-plant  population  density  studies,  poses  a  major  challenge. 
Cultivars  may  differ  greatly  in  their  yield  dynamics,  and  as  such  they  may  respond 
differently  to  changes  in  plant  population.   If  this  is  the  case,  then  the  analysis  of 
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combined  data  sets  from  different  cultivars  may  yield  results  which  misrepresent 
biological  and  ecological  realities. 


CHAPTER  5 
RESULTS  AND  DISCUSSION-1988  EXPERIMENT 


The  first  section  of  this  chapter  focuses  on  i)  the  effect  of  plant  population 
density  on  crop  variables  ii),  the  extent  of  the  correlation  of  these  variables  with 
seed  dry  weight  per  unit  area  and  iii),  the  yield  implications  of  i)  and  ii).  In  the 
second  section  the  cross  referencing  of  plant  population  density  effects  with  the 
strength  of  the  yield  association  of  variable,  is  used  to  establish  yield  determining 
relationships  (between  variables)  which  are  of  importance  in  yield-plant 
population  density  studies. 


Plant  Population  Density  Effects.  Yield  Correlations  and  their  Implications 

for  Yield 


Seed  dry  weight  (g  m"^)  increased  with  plant  population  density 
(Fig.  5-1).  Attempts  to  explain  this  result  were  done  by  categorizing  the  variables 
which  have  been  measured  as  being  i)  vegetative  ii),  related  to  nutrient  status 
and  iii),  reproductive.  The  effect  of  plant  population  density  on  R3,  R5  and  R8 
variables  was  determined  by  conducting  linear  regressions  on  data  which  were 
linearized  through  logarithmic  transformations.  These  models  only  show  if  plant 
population  density  had  a  significant  effect  on  the  variables  in  question,  and  with 
the  aid  of  the  model's  slope,  whether  this  effect  was  positive  or  negative. 
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F  ig .  5- 1 .  The  effect  of  p  la  nt  pop  u  latio  n  de  ns  ity  o  n  seed  a  nd  pod 
weight  of  Kirby  at  maturity,  1 988. 
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Correlations  of  variables  with  seed  dry  weight  (g  m  '^)  were  determined  for 
untransformed  data. 

Vegetative  Structures 

Drv  weight  response  to  plant  population  density 

At  R3,  total  dry  weight  (g  m"^)  ,  leaf  dry  weight  (g  m'^)  and  stem  dry  matter 
(g  fn"^)  (Fig  5-2)  (Table  5-1),  were  all  positively  affected  by  plant  population 
density.  At  R5,  total  dry  weight  (g  m"^),  leaf  dry  weight  (g  m"^)  and  stem  dry 
weight  (g  m"^)  (Fig.  5-3)  (Table  5-2)  were  all  positively  affected  by  plant 
population  density.  These  results  are  consistent  with  those  of  other  researchers 
(Parvez  et  al.,  1989;  Blumenthal  et  al.,  1988;  Udoguchi  and  McCloud,  1986). 

At  R3  and  R5  total  dry  weight  per  plant  decreased  (Fig.  5-4)  and  was 
significantly  affected  by  plant  population  density  (Tables  5-1  and  5-2).  Total  dry 
weight  decreased  sharply  up  to  a  density  of  25  plants  m'^,  and  tended  to  level  off 
after  44.4  plants  m'2.  for  both  R3  and  R5  (Fig.  5-4).  Obviously  at  higher  densities, 
increases  in  plant  population  density  were  not  as  effective  in  reducing  plant  dry 
weight. 

Drv  weights  fg  m'^)  and  correlations  to  seed  vield. 

Total  dry  weight  (g  m"^),  leaf  dry  weight  (g  m"^)  and  stem  dry  weight 
(g  m"2)  at  R3  were  positively  and  highly  correlated  with  seed  dry  weight  (g  m'^) 
at  R8  (Table  5-3).  Similar  results  were  obtained  at  R5  where  total  dry  weight 
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Table  5-1 .  Statistical  models  for  soybean  plant  variables  (Y)  with  plant 
density  (X)  at  R3  stage,  1988. 

Variable  (Y)  Statistical  model  r2 

Total  dry  weight(g  m"^)  InY  =  0.41  InX  +  5. 1 4  0.86" 

Leaf  dry  weight  (g  n\-\  InY  =  0.34lnX  +  4.43  0.78" 

Stem  dry  weight  (g  m"^j  InY  =  0.47lnX  +  4.47  0.88" 

Total  dry  weight  plant  "'•(g)  InY  =  0.58lnX +  5.14  0.92" 

*  ,**  significance  at  the  5  and  1%  levels,  respectively. 


Table  5-2.  Statistical  models  for  soybean  plant  variables  (Y)  at  R5  stage  with 
plant  density  (X),  1988. 

Variable  (Y)  Statistical  model  r2 

Total  dry  weight  (g  m"2)  InY  =  0.43lnX  +  5.34  0.79" 

Leaf  dry  weight  (g  rn^  InY  =  0.41  InX  +  4.49  0.75" 

Stem  dry  weight  (g  m"^j  InY  =  0.41  InX +  4.78  0.56" 

Total  dry  weight  plant '■'(g)  InY  =  0.52lnX  +  5.28  0.87" 


',**  significance  at  the  5  and  1%  levels,  respectively. 
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Fig.  5-2.  The  effect  of  plant  population  density  on  Stem,  leaf  and 
total  dry  weights  of  Kirby  at  the  R3  stagej  988. 
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Fig.  5-3.  The  effect  of  plant  population  density  on  Stem,  leaf  and 
total  dry  weights  of  Kirby  at  the  R5  stage,  1988. 
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Fig.  5-4.  The  effect  of  plant  population  density  on  vegetative  dry 
weight  of  Kirby  at  the  R3  and  R5  stages,  1 988. 
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Table  5-3.  Correlationxoefficients  of  soybean  plant  variables  at  R3  stage  with 
seed  dry  weight  (g  m"^)  at  R8, 1988. 


Variable 


Total  dry  weight  (g  m"^)  0.80" 

Leaf  dry  weight  (g  m"^)  0.77" 

Stem  dry  weight  (g  m"^)  0.81 " 


*,  "  significance  at  the  5  and  1%  levels,  respectively. 


Table  5-4.  Correla;tjon  coefficients  of  soybean  plant  variables  at  R5  stage  with 
seed  weight  (g  m"^)  at  R8,  1988. 


Variable  (Y) 


Total  dry  weight  (g  m"^)  0.63" 

Leaf  dry  weight  (g  m'^j  0.65" 

Stem  dry  weight  (g  m'^)  0.54" 


',**  significance  at  the  5  and  1%  levels,  respectively. 
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(g  m*2),  leaf  dry  weight  (g  m'2)  and  stem  dry  weight  (g  m"2)  at  R5,  were 
moderately  correlated  with  seed  dry  weight  (g  m"^)  at  R8  (Table  5-4).  Udoguchi 
and  McCloud  (1986)  and  Parvez  et  al.  (1989)  also  had  similar  findings. 
However,  Basuchaudhari  (1987)  found  only  moderate  correlations. 


Yield  implications. 

Total  vegetative  dry  weight  per  plant  decreased  with  plant  population 
density  and  therefore,  there  was  less  potential  for  higher  yields  on  a  per  plant 
basis.  However,  the  increase  in  plant  population  offset  the  reduction  in  total  dry 
weight  per  plant,  and  as  a  result  total  dry  weights  (g  m"^)  increased  with  plant 
population  density.  The  high  correlation  of  dry  weight  (g  m"^)  with  seed  dry 
weight  (g  m  '^),  indicates  that  accompanying  such  dry  weight  increases  per  unit 
area  are  i)  increased  photosynthetic  capacity  per  unit  land  area  and  ii),  increased 
amounts  of  nutrients  per  unit  land  area  for  remobiiization.  A  better  appreciation 
of  the  role  of  dry  weight  in  productive  processes  comes  with  examining  measures 
of  its  distribution-specifically  leaf  fractions,  the  leaf  area  index,  etc. 

Leaf  fraction,  leaf  area  and  light  intercepted  and  linear  reoressions 

At  R3,  the  stem  fraction  was  positively  affected  while  leaf  fraction  was 
negatively  affected  by  plant  population  density  (Fig.  5-5)  (Table  5-5).  However, 
at  R5  plant  population  density  had  no  significant  effect  on  the  stem  and  leaf 
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Fig.  5-5.  The  effect  of  plant  population  density  on  stenn  and  leaf 
fractions  of  Kirby  at  the  R3  and  R5  stages,  1  988. 
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Table  5-5.  Statistical  models  for  soybean  plant  variables  (Y)  with  plant 
density  (X)  at  R3  stage,  1988. 


Variable  (Y) 


Statistical  model 


r2 


Stem  fraction 

Leaf  fraction 

Percent  ligtit  intercepted 

Leaf  area  index 

Specific  leaf  area  (m^  g"^) 

Specific  leaf  weight  (g  m  "^) 


InY  =  O.OSlnX  -  .66 

0.61" 

InY  =  -0.07lnX  -  .71 

0.63" 

InY  =  0.09lnX  +  4.29 

0.59" 

InY  =  0.40lnX  +  .63 

0.83" 

InY  =  0.07lnX  -  3.81 

0.25* 

InY  =  -0.07lnX  +  3.81 

0.25* 

significance  at  the  5  and  1%  levels,  respectively. 


Table  5-6  .  Statistical  models  for  soybean  plant  variables  (Y)  with  plant  density 
(X)  at  R5  stage,  1988. 


Variable  (Y) 


Statistical  model 


r2 


Percent  light  intercepted 

Stenn  fraction 

Leaf  fractbn 

Leaf  area  index 

Specific  leaf  area  (m^  g"  M 

Specific  leaf  weight(g  m"^) 


lnY  = 

0.07lnX  +  4.35 

0.56* 

lnY  = 

-0.02lnX  -  .56 

0.01 

lnY  = 

-0.02lnX  -  .85 

0.02 

lnY  = 

0.38lnX  +  .80 

0.67* 

lnY  = 

-0.03lnX  -  3.69 

0.01 

lnY  = 

0.03lnX  +  3.69 

0.01 

significance  at  the  5  and  1%  levels,  respectively. 
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fractions  (Fig.  5-5)  (Table  5-6),  although  the  plant  density  trends  were  similar  to 
those  at  R3.  At  R3,  specific  leaf  area  (m^  g-^)  (Fig.  5-6)  and  specific  leaf  weight 
(g  m"2)  were  affected  by  plant  population  density  (Table  5-6).  At  R5,  specific  leaf 
area  (m^  g  ^l)  (Fig.  5-6)  and  specific  leaf  weight  (g  m*^)  were  not  affected  by 
plant  population  density  (Table  5-6).  By  comparison,  Nelson  and  Schweitzer 
(1988),  observed  that  specific  leaf  weight  (g  m"^)  (and  hence,  specific  leaf  area 
(m^  g""! ) )  was  not  greatly  affected  by  environments. 

At  higher  plant  population  densities,  due  to  soybean's  phyllotaxy,  it  is 
expected  that  a  greater  proportion  of  leaves  will  be  shaded.  However,  the  lack  of 
differences  in  specific  leaf  weight  at  R5,  may  have  been  due  to  the  fact  that  there 
was  probably  greater  differences  within  canopies,  than  between  canopies.  The 
specific  leaf  weights  and  specific  leaf  areas  taken  were  average  measures  over 
the  whole  canopy  and  they  do  not  indicate  differences  with  height  in  canopies. 
At  R3,  the  leaf  area  index  (Fig.  5-7)  and  percent  light  intercepted  (Fig.  5-8),  were 
affected  by  plant  population  density  (Table  5-5).  At  R5,  both  the  leaf  area  index 
(  Fig.  5-7  ,Table  5-6  )  and  the  percent  light  intercepted  (  Fig.  5-8,  Table  5-6  )  were 
affected  by  plant  population  density.  Light  interception  reached  a  maximum  of 
99%  and  leveled  off  at  a  population  density  of  25  plants  m"^. 

Leaf  fraction,  leaf  area  and  light  Intercepted  and  correlations  with  seed  yield. 

At  R3,  stem  and  leaf  fractions  were  highly  correlated  with  seed  dry  weight 
(g  m'2)  at  R8  (Table  5-7).  However,  at  R5  stem  and  leaf  fractions  were  not 
correlated  with  seed  dry  weight  (g  m"^)  (Table  5-8)  .  At  R3,  specific  leaf  weight 
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Fig.  5-6.  The  effect  of  plant  population  density  on  specific  leaf  area 
of  Kirby  at  the  R3  and  R5  stages,  1  988. 
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Fig.  5-7.  The  effect  of  plant  population  density  on  the  leaf 
area  index  of  Kirby  at  the  R3  and  R5  stages,  1  988. 
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Fig.  5-8.  The  elTectof  plant  population  density  on  percent 
light  interception  of  Kirby  at  the  R3  and  R5  stages,  1 988. 
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Table  5-7.  Correlati 
with  seed  wei 


rrelation  coefficients  of  soybean  plant  variables  at  R3  stage 
ght(g  m-2)at  R8,  1988. 


Variable  r 

Stem  fraction  0.69" 

Leaf  fraction  -0.69" 

Percent  light  intercepted  0.67" 

Leaf  area  index  0.79" 

Specific  leaf  area  (m^g"^)  0.41* 

SJjecific  leaf  weight  (g  m  '^)  -0.46" 


*,**  significance  at  the  5  and  1%  levels,  respectively. 


Table  5-8.  Correlai'on  coefficients  of  soybean  plant  variables  at  R5  stage  with 
seed  weight  (g  m"^)  at  R8, 1988. 

Variable  r 

Percent  light  intercepted  0.65" 

Stem  tract bn  -0.11 

Leaf  fraction  0.07 

Leaf  area  index  0.65 

Specific  leaf  area  (m^  g"^l  -0.10 

Specific  leaf  weight  (g  m"^)  -0.10 


significance  at  the  5  and  1%  levels,  respectively. 
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(g  m"2)  and  specific  leaf  area  (m^  g  ■"•)  were  moderately  correlated  with  seed 
weight  (g  m'^)  at  R8.   At  R5  specific  leaf  weight(g  m'2)  and  specific  leaf  area 
(m^  g  -I)  were  poorly  correlated  with  seed  weight  (g  m'^)  at  RB. 

Leaf  area  index  and  percent  light  intercepted  at  R3,  were  highly 
correlated  with  seed  weight  (g  m"2)  at  R8  (Table  5-7).  Leaf  area  index  and 
percent  light  intercepted  at  R5  were  highly  correlated  with  seed  weight  (g  m"^)  at 
R8  (Table  5-8).  This  is  consistent  with  the  findings  of  Basuchaudhari  (1987). 
However,  Ramseur  et  al.  (1985)  and  Costa  et  al.  (1980)  found  that  LAI  was  not 
correlated  with  yield. 

Yield  implications 

At  R3  with  increasing  plant  population  density  the  stem  fraction  increased. 
However,  at  R5  the  plant  population  effect  on  the  stem  fraction  was  not 
significant.  Specific  leaf  weights  tended  to  be  poorly  correlated  with  seed  dry 
weight. 

Bunce  (1988)  indicated  that  even  though  LAI  increased  with  plant 
population  density,  the  photosynthetic  capacity  per  unit  leaf  area  decreased, 
because  of  the  reduced  photosynthetic  capacity  of  shaded  leaves.  The  lower 
respiration  rate  of  shaded  leaves  results  in  them  not  being  parasitic  (Shibles  and 
Weber,  1965;  Duncan  et  al.  1967).  However,  it  is  likely  that  the  increase  in  leaf 
photosynthetic  area  per  unit  land  area,  more  than  makes  up  for  reduced 
photosynthetic  capacity  per  unit  leaf  area,  thereby,  increasing  the  potential  for 
increased  yields  with  increased  plant  population  density. 
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Nutrient  Status  and  Linear  Regression. 

At  R5,  plant  population  density  had  no  effect  on  N  concentrations  in  the 
stems  and  leaves  (Fig.  5-9  and  Table  5-9  ).  However,  N  content  (g  m"2)  in  the 
stems  and  leaves  were  affected  significantly  (Fig.  5-10,  Table  5-9  ).  At  R5,  plant 
population  density  had  no  effect  on  P  concentrations  in  either  stems  or  leaves  ( 
Fig.  5-11  and  Table  5-9  ).  Phosphorus  content  (g  m'^)  in  the  stems  and  leaves 
were  affected  significantly  (  Fig.  5-12  and  Table  5-9).  These  results  corroborate 
those  of  Abrahamson  and  Caswell  (1982)  and  Hanway  and  Weber  (1971). 
Unlike  the  findings  of  Buttery  (1969),  N  concentration  did  not  decrease  because 
of  inter-plant  competition  for  limited  N  resources.  Nitrogen  (g)  per  unit  leaf 
area  (  Fig.  5-13,  Table  5-9  ),  and  P  (g)  per  unit  leaf  area  (  Fig.  5-14,  Table  5-9) 
were  not  affected  by  plant  population  density.  Leaf  N  (g  m"^)  appeared  to  be 
greater  in  the  leaves  than  in  the  stems,  (the  intercept  and  slopes  of  the  statistical 
models  for  N  (g  m"^)  and  P  (g  m"^)  imply  this)  and  this  observation  is  consistent 
with  the  findings  of  Wang  and  Liu  (1987). 

Nutrient  status  and  correlations. 

The  correlation  results  are  provided  in  Table  5-10.  The  concentrations  of 
N  in  both  stems  and  leaves  at  R5  were  not  correlated  with  seed  dry  weight 
(g  m'2)  at  R8.   The  concentration  of  P   in  both  stems  and  leaves  at  R5  were  not 
correlated  with  seed  dry  weight  (g  m'^)  at  R8.  N  and  P  content  per  unit  leaf  area 
was  not  correlated  with  seed  dry  weight,  and  therefore  were  not  a  basis  for 
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Fig.  5-9.  The  effect  of  plant  population  density  on  leafandstenn 
N  concentrations  of  Kirby  at  the  R5  stage,  1 988. 
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Table  5-9  .  Statistical  models  for  soybean  plant  variables  (Y)  at  R5  stage  with 
plant  density  (X),  1988. 


Variable  (Y) 


Statistical  model 


r2 


1, 


Leaf  N  concentration 

Stem  N  concentration  (g  kg" ') 

Leaf  P  concentration  (g  kg'^J 

Stem  P  concentration  (g  kg' ' 

Leaf  N  (g  m'^) 

Leaf  N  (g  m'^  (leaf  area )) 

Stem  N  (g  m"^) 

Leaf  P  (g  m'^) 

Leaf  P  (g  m"^  (leaf  area ) 

Stem  P  (g  m"^) 


InY  =  0.04lnX  -  3.35  0.03 

InY  =  -0.02lnX  -  4.00  0.01 

InY  =  -0.03lnX  -  5.55  0.02 

InY  = -0.01  InX- 5.71  0.01 

InY  =  0.45inX  +  5.74  0.72" 

InY  =  0.07lnX  +  4.95  0.07 

InY  =  0.39lnX  +  5.38  0.45" 

InY  =0  .38lnX  +  3.55  0.64" 

InY  =  -0.01  InX +  2.75  0.00 

InY  =  0.40lnX  +  3.68  0.49" 


significance  at  the  5  and  1%  levels,  respectively. 


Table  5-10.  Correlatioricoefficients  of  plant  variables  at  R5  stage 
with  seed  weight  (g  m"^)  at  R8, 1988. 


Variable 


Leaf  N  concentration  (g  kg"^J 

Stem  N  concentration  (g  kg   ) 

Leaf  P  concentration  (g  kg"^ ) 

Stem  P  concentration  (g  kg"') 

Leaf  N  (g  m"^) 

Leaf  N  (g  m"^  (leaf  area )) 

Stem  N  (g  m"^) 

Leaf  P  (g  m'^) 

Leaf  P  (g  m'^  (leaf  area ) 

Stem  P  (g  m"^) 


0.31 
-0.02 
-0.13 
-0.05 
0.70" 
0.30 
0.51" 
0.64" 
0.02 
0.52" 


significance  at  the  5  and  1%  levels,  respectively. 
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Fig.  5-10.  The  effect  of  plant  population  density  on  leaf  and  stem 
N  content  of  KIrby  at  the  R5  stage,  1 988. 
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Fig.  5-11.  The  effect  of  plant  population  density  on  leaf  and  stem  P 
concentrations  of  KIrby  at  the  R5  stage,  1 988. 
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Fig.  5-12.  The  effect  of  plant  population  density  on  leaf  and  stem  P 
content  of  Kirby  at  the  R5  stage,  1  988. 
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Fig.  5-13.  The  effectofplantpopulation  density  on  N  content  per 
unitareaof  Kirbyatthe  R5  stage,  1988. 
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Fig.  5-14.  The  effect  of  plant  population  density  on  Pcontentper 
unit  leaf  area  of  Kirby  at  the  R5  stage,  1988. 
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differing  photosynthetic  capacity.  However,  leaf  N  and  P  content  were  highly 
correlated  with  seed  dry  weight  (g  m"^)  at  R8,  while  stem  N  and  P  content  were 
moderately  correlated  with  seed  dry  weight  (g  m'^)  at  R8. 

Yield  implications. 

There  were  no  apparent  deficiencies  of  N  and  P  (Figs.  5-9  and  5-11).  The 
N  concentration  was  comparable  to  values  reported  for  well  nodulated  crops.  It 
has  been  shown  that  P  and  N  concentrations  in  leaves  are  related  to 
photosynthetic  capacity  (Boote  et  al.,  1978;  Evans, 1983;  Sawada  et  al.,  1983; 
Terry  and  Ulrich,  1973)  and  hence  yield.  The  fact  that  plant  population  density 
had  no  effect  on  N  and  P  concentrations  in  both  leaves  and  stems  does  not 
diminish  the  importance  of  N  and  P  concentration  in  yield  determination. 

Total  N  and  P  uptake  on  the  other  hand,  like  yield,  are  largely  under 
environmental  control,  because  each  is  expressed  on  a  per  unit  land  area  basis 
in  a  plant  population  density  study.  Nitrogen  and  P  content  results  indicate  i),  a 
greater  remobilization  potential  at  higher  plant  population  densities,  as  long  as 
the  relative  rate  of  nutrient  translocation  to  seeds  remains  constant  and  ii),  a 
larger  photosynthetic  machinery  per  unit  land  area  at  higher  population 
densities.  These  two  events  can  explain  increased  yield  with  increasing  plant 
population  density. 
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Table  5-1 1 .  Statistical  models  for  soybean  plant  variables  (Y)  at  R5  stage  with 
plant  density  (X),  1988. 


Variable  (Y) 


Statistical  model 


r2 


Pod  dry  weight  (g  m''^) 
Pod  number  m"^ 


InY  =  0.52lnX  +  1 .76  0.61" 

InY  =  0.63lnX  +  4.79  0.71" 


',**  significance  at  the  5  and  1%  levels,  respectively. 


Table  5-12.  Statistical  models  of  plant  density  (X)  with  soybean  yield 
components  (Y)  at  R8,  1988. 


Variable  (Y) 


Statistical  model 


R^ 


Seed  dry  weight  (g  m"2) 
Pod  dry  weight  (g  m"2) 
Seed  number  m'2 
Pod  number  m"2 
Seed  number  Pod  "1 
Weight  seed "''  (g) 
Seed  dry  weight  plant  "1  (g) 
Pod  dry  weight  plant  '"^  (g) 
Seed  number  plant  "1 
Pod  number  plant  ""I 


InY  =  0.47lnX  +  3.73 
InY  =  0.46lnX  + 4.14 
InY  =  0.51lnX  + 6.21 
InY  =  0.49lnX  +  5.53 
InY  =  0.02lnX  +  .68 
InY  =  -0.04lnX  -  2.49 
InY  =  -0.53lnX  +  3.73 
InY  =  -0.54lnX  + 4.14 
InY  =  -0.49lnX  +  6.22 
InY  = -0.51  lnX  + 5.53 


0.75" 

0.80" 

0.93" 

0.91" 

0.05 

0.04 

0.79" 

0.84" 

0.92" 

0.92" 


significance  at  the  5  and  1%  levels,  respectively. 
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Reproductive  Structures 

Linear  regression  of  reproductive  structures  with  plant  population  density 

Pod  dry  weight  (g  m'2)  and  pod  number  (  m"2)  at  R5  were  significantly 
affected  by  plant  population  density  (Table  5-11).  Therefore,  it  was  not  surprising 
that  these  two  variables  were  also  shown  to  be  affected  by  plant  population 
density  at  R8;  pod  dry  weight  (g  m'^)  (Fig.  5-1,  Table  5-12)  and  pod  number 
(  m"2)  (Fig.  5-15,  Table  5-12)  Seed  dry  weight  (g  m"2)  was  affected  positively 
by  plant  population  density  (  Fig.  5-1,  Table  5-12).  This  result  is  consistent  with 
the  results  of  Parvez  et  al.  (1989),  Spilde  et  al.  (1980),  Doss  and  Thurlow  (1974) 
and  Caviness  (1966). 

Seed  number  m"^  at  R8  increased  with  plant  population  density  (  Fig. 
5-15,  Table  5-12).  This  result  corroborates  the  findings  of  Leuschen  and  Hicks 
(1977)  and  Pontes  and  Ohirogge  (1972).  Weight  seed  "''(g)  was  not  affected  by 
plant  population  density  (Fig  5-16).  Seed  number  pod"""  was  not  affected  by 
plant  population  density  (Fig.  5-17).  Dominguez  and  Hume  (1978),  and  Pontes 
and  Ohirogge  (1972)  had  similar  results.  However,  this  finding  was  contrary  to 
those  of  Johnson  and  Harris  (1967)  (with  the  soybean  cultivar  'Lee'),  Parks  and 
Manning  (1980),  Wilcox  (1974),  Wright  et  al.  (1987)  and  Weil  and  Ohirogge 
(1976). 

Seed  weight  and  pod  dry  weight  per  plant  (Fig.  5-18)  decreased  with 
plant  population  density,  and  were  also  significantly  affected  (Table  5-12).   Pod 
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Fig.  5-15.  The  effect  of  plant  population  density  on  seed  and  pod 
number  of  Kirby  at  maturity,  1 988. 
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Fig.  5-16.  The  effectof  plant  population  density  on  vy^ight  per 
seed  of  Kirby  at  maturity,  1988. 
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Fig.  5-1 7.  The  effect  of  plant  population  on  seed  number  per 
pod  of  Kirby  at  maturity,  1 988. 
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Fig.  5-1 8.  The  effect  of  plant  population  density  on  seed  and 
pod  weight  per  plant  of  Kirby  at  maturity,  1988. 
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and  seed  number  per  plant  responded  in  a  similar  manner  (Fig.  5-19,  Table 
5-12).  All  variables  tending  to  level  off  at  above  44.4  plants  m'^. 

Correlations  of  reproductive  structures  with  vield 

Pod  dry  weight  (g  m"^)  and  pod  number  (m"2)  at  R5  were  both  highly 
correlated  with  seed  dry  weight  (g  m'^)  at  R8  (Table  5-13).  Among  R8  variables 
(Table  5-14),  pod  dry  weight  (g  m"^),  seed  number  (m"^)  and  pod  number 
(m"2)  were  all  highly  correlated  with  seed  dry  weight  (g  m"^)  (Table  5-14).  These 
results  are  consistent  with  those  of  Basuchaudhari  (1987),  Kaw  and  Menon 
(1972),  Pontes  and  Ohirogge  (1972),  Parvez  et  al.  (1989).  Seed  number  per  pod 
and  weight  per  seed  (g)  were  not  correlated  with  seed  dry  weight  (g  m'^)  (Table 
5-14).  These  findings  are  consistent  with  those  of  Board  (1985),  Tandojam 
(1986),  Johnson  et  al.  (1955),  and  Hartwig  and  Edwards  (1970).  However  Ma 
(1946)  found  that  seed  per  pod  was  correlated  with  seed  dry  weight  (g  m'^),  and 
Weber  and  Moorthy  (1952)  reported  that  there  was  a  high  positive  correlation 
between  weight  per  seed  and  seed  dry  weight  (g  m"2)  (Table  5-14). 

Yield  implications 

The  reproductive  variables  which  can  account  for  increased  yield  are  pod 
and  seed  number  m  '^.  Variables  unable  to  account  for  yield  increases  are 
weight  seed  '^  and  seed  pod  "^.  Unlike  the  former  two  variables,  the  latter  two 
variables,  because  they  are  not  expressed  on  a  per  area  basis  are  not  easily 
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Fig.  5-1 9.  The  effect  of  plant  population  density  on  seed  and 
pod  number  per  plant  of  Kirby  at  maturity,  1988. 
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Table  5 
with  seed 


13.  Correlation  coefficients  of  soybean  plant  variables  at  R5  stage 
jd  weight  (g  m"^)  at  R8, 1988. 


Variable  (Y) 


Pod  dry  weight  (q  m"2)  0.61 ' 

Pod  number  m"^  0.71' 


*,**  significance  at  the  5  and  1%  levels,  respectively. 


Table  5-14.  Correlation  coefficients  of  soybean  plant  variables  at  R8  with 

seed  weight  (g  m*^)  at  R8, 1988. 

Variable  (Y)  r 

Pod  dry  weight  (g  nr^)  0.99" 

Seed  number  (m-2)  0.92" 

Pod  number  m-2  0.90" 

Seed  number  pod  ""^  0.32 

Weight  seed  -''  (g)  0.23 


*,**  significance  at  the  5  and  1%  levels,  respectively. 
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subject  to  changes  in  plant  population  density.  It  should  be  emphasized  that 
weight  seed  "^(g),  and  seed  pod  '^ ,  make  vital  contributions  to  final  yield,  but 
they  are  unable  to  account  for  the  plant-density-mediated  variations  in  yield. 

Results  and  Discussions  of  Cross  Tables 

In  this  sub-section,  the  relationship  between  variables  will  be  examined 
using  the  following  'degree-of-correlation'  classification  system: 

Highly  correlated  with  yield  or  affected  by  plant  population  density:  r  or  R^  : 
±0.75-  1.00.  (H) 

Moderately  correlated  with  yield  or  affected  by  plant  population  density:  r  or  R^  : 
±  0.50  -  0.74  (M) 

Low  correlation  with  yield  or  not  affected  by  plant  population  density:  r  or  R^  : 
±  0.00  -  0.49  (L) 

Yield  and  its  components  at  R8 

The  variables  at  R8  which  are  highly  correlated  with  seed  dry  weight 
(g  m"2),  are  the  same  variables  which  are  greatly  affected  by  plant  population 
density  (Table  5-15).  Among  the  variables  highly  associated  with  yield  and 
greatly  affected  by  plant  population  density,  the  relationship  between  pod 
number  (m"2)  and  seed  number  (m"^)  appears  to  be  the  strongest.  This 
observation  can  be  attributed  to  two  facts.   First,  both  variables  are  expressed  on 
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Table  5-15.  Yield-Plant  Density  Cross-Table  for  soybean  plant  variables  highly 
associated  with  yield  and  exhibiting  high  plant  density  effect,  at  R8,  1988. 


Variable 


High  yield  association  (R8) 


High  plant  density  effect  (R8) 


Pod  dry  weight  (g  m"2) 


Seed  number  (m'2) 
Pod  number  (m'2) 


0.94" 
0.92" 


Seed  number  (m"2) 
Pod  number  (m'2) 


Pod  number  (m"2) 
Pod  dry  weight  (g  m"2) 

Seed  number  (m'2) 
Pod  dry  weight  (g  m"2) 


0.94" 
0.94" 

0.98" 
0.92" 


*,"  significance  at  the  5  and  1%  levels,  respectively. 


104 


a  land  area  basis.  Secondly,  pod  number  and  seed  number  are  related  via 
seeds  pod  "^,  and  it  has  been  shown  that  seed  pod"^  is  not  affected  by  plant 
population  density.  Therefore  seed  number  and  pod  number  will  tend  to  vary 
with  each  other,  resulting  in  a  high  r  value.  All  other  relationships  between  pod 
dry  weight  (g  m"^),  pod  number  (m"2)  and  seed  number  (m*2)  also  seem  to  be 
the  primary  relationships  of  great  importance  in  yield  determination.  Again,  this 
observation  in  part,  can  be  attributed  to  the  fact  that  these  variables  are 
expressed  on  a  per  area  basis. 

The  poor  relationships  between  yield  components  not  affected  by  plant 
population  density,  namely  seed  pod'^  and  weight  seed'^,  and  components 
highly  associated  with  yield  (Table  5-16),  reinforces  the  point  that  seed  pod"'' 
and  weight  seed'"'  are  far  removed  from  mediating  plant  population  density 
effects  on  yield.  It  is  interesting  to  note  that  seed  pod'^,  and  weight  seed"^,  are 
poorly  correlated  (Table  5-17).  Therefore  even  if  there  are  more  seeds  pod'^, 
weight  seed"""  tends  to  remain  constant. 

At  R5,  the  variables  total  dry  weight  (g  m"^)  and  leaf  dry  weight  (g  m"^), 
which  exhibit  high  plant  population  density  effect  are  moderately  correlated  with 
pod  dry  weight  (g  m'2),  seed  number  (m"^)  and  pod  number  (m"2)  at  R8  (Table 
5-18).  The  physiological  basis  of  this  is  to  be  found  in  i),  greater  vegetative  dry 
weight  results  in  a  greater  capacity  to  remobilize  nutrients,  and  ii),  leaf  dry  weight 
is  a  reflection  of  the  photosynthetic  capacity  of  a  crop. 

Table  5-19  lists  variables  which  are  not  associated  with  yield  and  are  not 
affected  by  plant  population  density.  All  the  variables  except  for  stem  N  and  P 
contents  are  expressed  in  units  other  than  on  a  land  area  basis,  and  this  may 
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Table  5-16.  Yield-Plant  Density  Cross-Table  for  soybean  plant  variables  highly 
associated  with  yield  and  exhibiting  no  plant  density  effect  at  R8, 1988. 


Variable 


High  yield  association  (R8) 


No  plant  density  effect  (R8) 


Pod  dry  weight  (g  m"2) 
Seed  number  (m"2) 
Pod  number  (m"2) 


Seed  pod'"" 
Weight  seed '"I  (g) 

Seed  pod"'' 
Weight  seed  ""'(g) 

Seed  pod'"" 
Weight  seed  ""'(g) 


0.31 
0.17 

0.25 
-0.14 

0.08 
-0.17 


*,**  significance  at  the  5  and  1%  levels,  respectively. 
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Table  5-17.  Yield-Plant  Density  Cross-Table  for  soybean  plant  variables  not 
associated  with  yield  and  exhibiting  no  plant  density  effect  at  R8, 1988. 


Variable 

No  yield  association  (R8) 

No  plant  density  effect  (R8) 

Seed  pod'"" 

Weight  seed  "'•(g) 

0.23 


*,**  significance  at  the  5  and  1%  levels,  respectively. 
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Table  5-18.  Yield-Plant  Density  Cross-Table  for  R5  soybean  plant  variables 
exhibiting  high  plant  density  effect  and  high  yield  association,  1988. 


High  yield  association 


Variable 


High  plant  density  effect 


Pod  dry  weight(R8)(g  m-2) 


Seed  number  (R8)(m"2) 


Pod  number{R8)(m"2) 


Total  dry  weight(R5)(g  m"2) 
Leaf  dry  weight(R5)(g  m"2) 

Total  dry  weight(R5)(g  m'2) 
Leaf  dry  weight(R5)(g  m"2) 

Total  dry  weight(R5)(g  m'2) 
Leaf  dry  weight(R5)(g  m'2) 


0.68" 
0.69** 

0.75** 
0.77** 

0.69** 
0.74** 


*,**  significance  at  the  5  and  1%  levels,  respectively. 
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Table  5-19.  Yield-Plant  Density  Cross-Table  for  soybean  plant  variables  not 
associated  with  yield  and  not  affected  by  plant  density  at  R5  and  R8, 1988. 

Variable 


No  yield  association 


No  plant  density  effect 


Leaf  N  concentration  (R5)  (g  kg'1] 
Stem  N  concentration  (R5)  (g  kg"') 
Leaf  P  concentration  (R5)  (g  kg""") 
Stem  P  concentration  (R5)  (g  kg"1) 
Stem  fraction  (R5) 
.  Leaf  fractbn  (R5) 
Specific  leaf  area  (R5)  (cm^  g-1) 
Specific  leaf  weighit  (R5)  (g  cm"2) 
Leaf  N  (R5)  (g  m-2  (leaf  area)) 
Leaf  P  (R5)  (g  m-2  (leaf  area)) 
Seed  number  Pod'"'  (R8) 
Weight  seed  ■■•  (R8)(g) 


Leaf  N  concentration  (R5)  (g  kg-"!) 

Stem  N  concentration  (R5)  (g  kg-"") 

Leaf  P  concentration  (R5)  (g  kg-1) 

Stem  P  concentration  (R5)  (g  kg-'') 

Stem  fraction  (R5) 

Leaf  fraction  (R5) 

Specific  leaf  area  (R5)  (cm^  g-^) 

Specific  leaf  weigfit  (R5)  (g  cm-2) 

Leaf  N  (R5)  (g  m-2  (leaf  area)) 

Stem  N  (R5)  (g  m-2) 

Leaf  P  (R5)  (g  m-2  (leaf  area)) 

Stem  P  (R5)  (g  m-2_| 

Seed  number  Pod-'(R8) 

Weight  seed -■'(R8)(g) 
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explain  some  of  the  unresponsiveness  of  these  variables.  Table  5-20,  provides 
some  interesting  relationships  between  variables  not  associated  with  yield  and 
not  affected  by  plant  population  density.  At  R5,  P  concentration  is  moderately 
and  positively  correlated  to  specific  leaf  area,  this  may  be  due  to  the  fact  that  P 
plays  an  important  in  leaf  expansion.  Also,  leaf  N  and  P  per  unit  leaf  area  are 
highly  and  positively  correlated  to  specific  leaf  weight.  Therefore,  thicker  leaves 
should  be  expected  to  have  a  larger  photosynthetic  apparatus. 

At  R3,  relationships  which  have  considerable  potential  for  mediating  plant 
population  density  effects  on  yield  are  total  dry  weight,  leaf  dry  weight,  stem  dry 
weight  and  leaf  area  index  (Table  5-21).  Since  specific  leaf  weight  is  unaffected 
by  plant  population  density,  it  should  not  be  surprising  that  leaf  area  (as  LAI)  is 
closely  associated  with  leaf  dry  weight.  The  remaining  relationships,  though 
important,  provide  nothing  revealing  as  most  are  directly  related,  such  as  leaf  dry 
weight  and  total  dry  weight. 

A  summary  and  the  relevance  of  cross-referencing  in  yield-plant 
population  density  studies  is  given  in  Table  5-22,  where  it  is  clearly  shown  that  all 
reproductive,  photosynthetic  and  remobilization  parameters  which  have  been 
found  to  be  plastic  are  associated  with  seed  weight  (g  m'^).  In  conclusion  the 
variables  which  through  their  relationships  have  the  greatest  potential  for 
explaining  plant-density-mediated  effects  on  yield  are  seed  number  (R8)  (m"^), 
pod  number  (R8)  (m"2),  pod  dry  weight  (R8)  (g  m"2),  total  dry  weight  (R5)  (g  m'2), 
and  leaf  dry  weight  (R5)  (g  m"^).  The  former  two  owe  their  importance  to  yield 
because  they  are  yield  components  which  are  expressed  on  land  area  bases; 
pod  dry  weight  because  it  incorporates  yield;  and  total  dry  weight  (R5)  and  leaf 
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Table  5-20.  Yield-Plant  Density  Cross-Table  for  soybean  plant  variables  not 
associated  with  yield  and  exhibiting  no  plant  density  effect  at  R5  and  which  show 
high  association  to  each  other,  1988. 


Variable 


No  yield  association 


No  plant  density  effect 


Leaf  N  concentration  (R5)  (g  kg""! 
Leaf  N  (R5)  (g  m-2  (leaf  area)) 
Leaf  N  (R5)  (g  m-2  (leaf  area)) 
Leaf  N  (R5)  (g  m-2  (leaf  area)) 
Leaf  P  concentration  (R5)  (g  kg'1 
Leaf  P  concentration  (R5)  (g  kg-"" 
Leaf  P  concentration  (R5)  (g  kg-1 
Stem  fraction  (R5) 
Specific  leaf  weiglit  (R5)  (g  m-2) 
Specific  leaf  weight  (R5)  (g  m-2) 
Specific  leaf  weight  (R5)  (g  m-2) 
Specific  leaf  weight  (R5)  (g  m-2) 
Leaf  fraction  (R5) 
Specific  leaf  area  (R5)  (m2  g-1) 
Specific  leaf  area  (R5)  (m2  g-1  j 
Specific  leaf  area  (R5)  (R5)  (m^  g 
Specific  leaf  area  (R5)  (m2  g-1) 
Stem  P  concentration  (R5)  (g  kg- 
Leaf  P  (R5)  (g  m-2  (leaf  area)) 
Leaf  P  (R5)  (g  m-2  (leaf  area)) 
Stem  N  concentration  (R5)  (g  kg- 
Stem  N  concentration  (R5)  (g  kg 
Stem  N  concentration  (R5)  (g  kg 


Stem  N  concentration  (R5)  (g  kg'"") 
Leaf  P  (R5)  (g  m-2  (leaf  area)) 
Specific  leaf  weight  (R5)  (g  m-2) 
Specific  leaf  area  (R5)  (m2  g-l) 
Specific  leaf  area  (R5)  (m2  g-lj" 
Specific  leaf  weight  (R5)  (g  m-^) 
Stem  N  concentration  (R5)  (g  kg-^) 
Leaf  fraction  (R5) 
Specific  leaf  area  (R5)  (m2  g-1) 
Leaf  P  (R5)  (g  m-2  (leaf  area)) 
Leaf  N  (R5)  (g  m-2  (leaf  area)) 
Leaf  P  concentration  (R5)  (g  kg-1) 
Stem  fraction  (R5) 
Specific  leaf  weight  (R5)  (g  m-2) 
Leaf  P  (R5)  (g  m-2  (leaf  area)) 
Leaf  N  (R5)  (g  m-2  (leaf  area)) ) 
Leaf  P  concentration  (R5)  (g  kg-1] 
Stem  N  concentration  (R5)  (g  kg-') 
Specific  leaf  weight  (R5)  {a  m-2) 
Specific  leaf  area  (R5)  (m^  g-1) 
Stem  P  concentration  (R5)  (g  kg-1) 
Leaf  N  concentration  (R5)  (g  kg-"") 
Leaf  P  concentration  (R5)  (g  kg-"") 


0.62" 

0.71" 

0.70" 

-0.64" 

0.60" 

0.54" 

0.57" 

-0.99" 

-0.93" 

0.78" 

0.70" 

-0.54" 

-0.99" 

-0.93" 

-0.69" 

-0.64" 

0.60" 

0.84" 

0.78" 

-0.69" 

0.84" 

0.62" 

0.57" 


significance  at  the  5  and  l%levels,  respectively. 
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Table  5-21 .  Yield-Plant  Density  Cross-Table  for  soybean  plant  variables  highly 
associated  with  yield  and  exhibiting  high  plant  density  effect  at  R3, 1988. 


Variable 


High  yield  association 


High  plant  density  effect 


Stem  dry  weight  (g  m"2) 


Total  dry  weight  (g  m'2) 


Leaf  area  index 


Leaf  dry  weight  (g  m"2) 


Stem  dry  weight  (g  m"2) 
Total  dry  weight  (g  m'2) 
Leaf  area  index 
Leaf  dry  weight  (g  m"^) 

Stem  dry  weight  (g  m"2) 
Total  dry  weight  (g  m"2) 
Leaf  area  index 
Leaf  dry  weight  (g  m'^) 

Stem  dry  weight  (g  m"2) 
Total  dry  weight  (g  m"^) 
Leaf  area  index 
Leaf  dry  weight  (g  m"2) 

Stem  dry  weight  (g  m"2) 
Total  dry  weight  (g  m"2) 
Leaf  area  index 
Leaf  dry  weight  (g  m"2) 


1 .00" 
0.99** 
0.97** 
0.98** 

0.99** 
1.00** 
0.97** 
0.99** 

0.97** 
0.97** 
1.00** 
0.96** 

0.98** 
0.99** 
0.96** 
1.00** 


*,  **  Significance  at  the  5  and  1%  levels,  respectively. 
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Table  5-22.  Plastic  and  non-plastic  crop  responses  of  reproductive, 
remobilization  and  photosynthetic  parameters  and  their  impact  on  seed  dry 
weight  (g  m"2),  1988. 

Reproductive  Parameters 


'V 


Plastic  responses 
Pod  number  (R5)(m"2) 
Pod  dry  weight(R5)(g  m' 
Pod  number  (RSlCm'^) 
Seed  number  (R8)(m"^) 
Pod  dry  weight(R8)(g  m" 


Non-plastic  responses 
Seed  number  Pod  "'(R8) 
Weight  seed '"'(RB) 


Remobilization  Parameters 


Plastic  responses 
Total  dry  weight  (R5)(g  m"^) 
Leaf  dry  weight  (R5)(g  m'^i 
Stem  dry  weight  (R5)(g  m"^) 
Leaf  N  (R5)(g  m'^) 
Leaf  P  (R5)(g  m'^) 
Stem  N  (R5)  (g  m'^) 
Stem  P  (R5)  (g  m'^) 


Non-plastic  responses 
Leaf  N  concentration  (R5)  (g  kg"') 
Leaf  P  concentration  (R5)  (g  kg'^| 
Stem  N  concentration  (R5)  (g  kg"') 
Stem  P  concentration  (R5)  (g  kg'^) 


Photosynthetic  Parameters 


Plastic  responses 
•V         Leaf  area  index  (R5) 

*  V  Leaf  dry  weight  (R5)  (g  m"^) 

*  V  Percent  light  intercepted(R5) 

*  V  Leaf  N  (R5)  (g  m'^) 

*  V  Leaf  P  (R5)  (g  m'^) 


Non-plastic  responses 
Specific  leaf  area  (R5)(m2  g""') 
Specific  leaf  weight(R5)(g  m'^) 
Leaf  N  concentration  (R5)  (g  kg"^) 
Leaf  P  concentration  (R5)  (g  kg"^) 
Leaf  N  (R5)  (g  m"2(ieaf  area) 
Leaf  fraction  (R5) 


*,  V,  responses  affected  by  plant  density  and  associated  with  seed  weight  (g  m"^), 
respectively. 
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dry  weight  (R5)  because  of  their  relevance  to  the  processes  of  photosynthesis 
and  remobilization. 


CHAPTER  6 
AN  INVESTIGATION  INTO  THE  YIELD-PLANT  DENSITY  RESPONSE  WITH 

SOYGRO 


In  this  chapter  the  process-oriented  crop  growth  simulator  for  soybean, 
SOYGRO,  is  used  to  investigate  Duncan's  following  assumptions  on  yield-plant 
density  relations: 

I)      that  percent  light  interception  at  growth  stage  R5     provides  an 
effective  basis  for  determining   and  comparing  the  seed  yield 
efficiencies  of  crops  of  different  plant  densities.  Seed  yield  efficiency 
is  defined  as  seed  yield  divided  by  percent  light  interception  at  R5. 

ii)  that  the  contribution  of  remobilization  to  seed  yield  has  no  differential 
effect  on  the  calculation  of  seed  yield  efficiencies  across  different 
planting  densities,  and 

Hi)  that  with  the  onset  of  interplant  competition,  there  is  an  increased 
efficiency  with  which  a  crop  uses  light  for  seed  production. 

Duncan's  theory  relies  in  large  part  on  these  three  assumptions.  First,  by 
relying  on  percent  light  interception  at  R5  as  a  basis  for  calculating  seed  yield 
efficiency,  Duncan  assumed  that  the  amount  of  light  intercepted  between  R5  and 
R7,  is  directly  proportional  to  the  percent  of  light  intercepted  at  R5.  For  example, 
in  Phase  II,  where  different  plant  densities  intercept  an  equivalent  95%  of  the 
incoming  light,  it  is  assumed  that  each  plant  density  intercepts  the  same  amount 
of  light  between  R5  and  R7.  This  was  the  basis  for  Duncan's  statement  that,  in 
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Phase  II,  if  yield  per  unit  area  increases  (and  since  percent  light  interception  is 
constant)  then  seed  yield  efficiency  (based  on  light  interception)  has  to  increase. 

The  findings  of  Shibles  and  Weber  (1965)  that  soybean  crop  growth  rate 
per  unit  ground  area  should  be  maximized  and  be  constant  when  the  canopy 
exhibits  complete  insolation  interception,  may  appear  to  support  Duncan's 
assumption.  However,  this  statement  is  not  to  be  taken  as  meaning  that  total 
canopy  photosynthesis  (growth  rate)  between  R5  and  R7  is  at  a  maximum  and  is 
constant  over  a  given  range  of  plant  densities.  Given  that  canopy  structure,  leaf 
distribution  and  degree  of  shading  are  all  greatly  influenced  by  plant  density,  it  is 
hardly  likely  that  plant  densities  corresponding  to  constant  yield  per  area  will 
consist  of  crops  with  similar  canopy  structures.  When  the  impact  of  nutrient 
remobilization  and  its  effect  upon  leaf  loss  and  canopy  structure  are  considered, 
it  is  hardly  likely  that  the  light  intercepted  between  growth  stages  R5  and  R7  will 
constant  over  a  range  of  plant  densities  which  provide  maximum  yield. 

Finally  Duncan  defined  'seed  yield  efficiency'  as  the  efficiency  with  which 
intercepted  PAR  is  used  in  seed  production  between  R5  and  R7.  Duncan, 
therefore,  confined  the  determination  of  'seed  yield  efficiency',  and  thus 
differences  in  yield  between  plant  densities  to  the  amount  of  light  intercepted 
during  R5  to  R7.  From  Duncan's  perspective,  the  process  of  nutrient 
remobilization  which  commences  at  R5  was  seen  strictly  as  a  physiological 
process  whose  potential  is  determined  by  the  vegetative  weight  of  the  plant  at  R5. 
Remobilization,  although  seen  by  Duncan  as  a  process  impacting  upon  yield, 
was  not  seen  as  a  process  which  differentially  impacts  upon  the  'seed  yield 
efficiency'  equation.  Duncan's  disregard  for  remobilization  is  justified  if  i),  the 
proportion  of  remobilized  nutrients  from  vegetative  organs  to   seeds  is  constant, 
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and  also  if,  ii)  the  proportion  of  remobilized  nutrients  in  the  seed  is  constant  over 
all  densities  investigated.  If  true,  then  differences  in  seed  yield  efficiency  could 
not  be  attributed  to  differences  in  the  rate  of  nutrient  remobilization;  instead,  yield 
differences  would  be  due  to  differences  in  the  efficiency  with  light  intercepted 
between  R5  and  R7  is  used.  If  this  is  the  case  then  Duncan's  arguments  on 
increased  seed  yield  efficiency  in  Phase  II  would  hold  although  the  calculated 
seed  yield  efficiencies  would  be  artificially  high.  If  Duncan's  assumption  does  not 
hold,  then  the  process  of  remobilization  confounds  any  calculation  of  seed  yield 
efficiency  which  relies  strictly  on  light  interception  at  or  after  R5. 

Duncan's  third  assumption  was  that  between  Phase  I  and  Phase  II  there 
is  a  range  of  plant  densities  within  which  there  is  an  increasing  degree  of  mutual 
shading  as  density  increases.  The  increase  in  yield  per  unit  area  observed 
within  this  range  of  plant  densities  is  assumed  to  be  due  both  to  the  increase  in 
light  intercepted  and  to  an  increase  in  the  efficiency  with  which  light  is  used. 
However,  Duncan  provided  no  basis  for  his  view  that  with  the  onset  of  shading 
(the  end  of  Phase  I)  there  is  an  increased  efficiency  with  which  PAR  is  used  in 
seed  production. 

A  Crop  Simulation  Perspective  on  Duncan's  Assumptions 

Duncan's  "thought  experiment"  and  the  assumptions  he  made  have  yet  to 
be  critically  investigated.  Egli's  (1988)  research  efforts  were  directed  specifically 
towards  Duncan's  yield-response  theory.  However,  Egli  did  not  question  the 
validity  of  the  assumptions  underlying  Duncan's  theory;  instead,  he  accepted 
Duncan's  postulates,  and  obtained  results  which  clearly  showed  the  existence  of 
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Duncan's  three  Phases.  While  actual  data  are  ideal  in  arriving  at  conclusions 
concerning  physiological  and  ecological  processes  operating  in  yield-plant 
density  relationships,  crop-simulation  modeling  can  serve  as  a  useful  tool  in  i) 
developing  theories  and  hypotheses  and  ii),  testing  the  validity  of  theories.  One 
such  model  which  holds  much  potential  for  investigating  Duncan's  theory,  and  for 
the  development  of  hypotheses  which  address  yield-plant  density  responses  in 
soybean  is  SOYGRO  (Jones  et  al.,  1989). 

SOYGRO  is  a  process-oriented  soybean  crop  growth  simulator  developed 
at  the  University  of  Florida  by  an  interdisciplinary  research  team.  The  simulator 
predicts  daily  growth,  partitioning,  soil-crop  water  balance  and  seed  yield 
response  to  daily  weather  (solar  radiation,  maximum  and  minimum  temperature 
and  rainfall).  The  simulator  is  sensitive  to  row  and  plant  spacing,  plant 
population,  planting  date  and  cultivar  selection,  and  irrigation  management 
options. 

A  version  of  SOYGRO  has  recently  been  adapted  to  include  a  hedgerow 
photosynthesis  -  light  interception  subroutine,  which  operates  on  a  hourly  step 
within  the  model  (Boote  et  al.,  1989).  This  version  predicts  hourly  (  and  daily  ) 
light  interception  and  canopy  photosynthesis  as  a  function  of  the  day  of  year, 
latitude,  solar  radiation,  temperature,  rainfall,  leaf  area  index,  canopy  height  and 
canopy  width.  The  latter  three  variables  are  predicted  internal  to  the  model. 
Additionally,  these  are  external  model  parameters  which  modify  the  increase  in 
height  and  width  as  a  function  of  vegetative  stage,  temperature, 
photosynthetically  active  radiation  (PAR),  and  water  deficit.  The  hedgerow 
version  of  SOYGRO  makes  it  suited  to  evaluate  response  of  light  interception  and 
yield  to  plant  population.    This  version  of  SOYGRO  is  still  in  its  developmental 
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stage,  and  since  the  model  has  only  been  tested  against  row  spacing  data,  it  is 
also  desirable  to  test  its  response  to  equidistant  plant  population  studies  of  Egli 
(1988). 

Egli's  data  set  covered  an  adequate  range  of  plant  densities;  however,  his 
data  set  had  insufficient  dry  matter  accumulation,  light  interception  and  nitrogen 
status  data  for  the  investigation  of  yield-plant  density  responses,  and  the 
examination  of  Duncan's  theory.  SOYGRO  was  used  to  generate  light 
interception,  dry  matter  accumulation,  N  status  data,  as  well  as  predict  yield 
response  to  plant  density.  Models  themselves,  however,  are  only  as  valid  as  the 
assumptions  under  which  they  operate.  It  is  from  this  perspective,  that  the 
assumptions  in  SOYGRO  will  be  discussed  relative  to  its  predictions  of  yield-plant 
density  responses. 

Simulation  Exercise 

The  crop  simulation  exercise  was  conducted  on  yield  data  from  Egli's 
(1988)  study  on  the  effect  of  plant  density  on  soybean  yield.  This  experiment  was 
conducted  in  1986,  near  Lexington  Kentucky  (38°N  Lat.)  with  the  soybean 
cultivar  Cumberland  [fvlG  III]  which  has  an  indeterminate  growth  habit.  The  soil 
was  a  Maury  silt  loam  (fine,  mixed,  mesic  Typic  Paleudalf)  and  planting  was  done 
in  "wagon  wheel"  design  which  resulted  in  an  approximate  square  planting 
arrangement  consisting  of  17  plant  populations.  Yield,  irrigation  management, 
soil  and  weather  data  pertinent  to  Egli's  experiment  were  collected  prior  to 
attempts  at  calibrating  SOYGRO. 
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The  yield  (Table  6-1  )  irrigation  schedule  (Table  6-2)  and  weather  data 
were  provided  by  Egli.  Soil  parameters  were  determined  using  results 
fromanalyses  of  the  United  States  Soil  Conservation  Service  (1983,  1964)  and 
these  included  percent  slope,  horizon  thickness,  and  for  each  horizon,  the  bulk 
density,  relative  root  presence,  organic  carbon  concentration,  pH,  percent  sand, 
percent  silt  and  percent  clay.  Using  a  procedure  available  as  part  of  the  IBSNAT 
(International  Benchmark  Sites  Network  for  Agrotechnology  Transfer)  crop 
growth  models  the  saturated  water  content,  the  drained  upper  limit  ,  the  lower 
limit  of  plant-extractable  water  and  the  saturated  hydraulic  conductivity  were 
derived  from  percent  sand,  percent  silt,  percent  clay  and  organic  carbon. 
Characteristics  of  the  derived  soil  which  represent  the  Maury  silt  loam  used  by 
Egli  are  given  in  Table  6-3. 

The  model  was  calibrated  by  altering  genetic  and  crop  growth 
parameters.  The  data  provided  by  Egli  for  model  calibration  included:  seed  yield 
per  unit  area,  seed  yield  per  plant,  seed  number  per  unit  area  and  weight  per 
seed.  In  addition  to  these  yield  variables,  SOYGRO  simulated  values  were 
output  for  percent  light  intercepted  over  the  growing  season;  percent  light 
intercepted  both  at  noon  and  daily  at  growth  stage  R5;  (beginning  seed 
anywhere  on  the  plant)  total  light  intercepted  during  the  vegetative  and 
reproductive  phases  and  seasonal  light  intercepted.  The  maximum  vegetative 
nitrogen  content,  mobilized  nitrogen  and  nitrogen  content  of  seeds  were  also 
predicted. 

No  previous  attempts  had  been  made  to  model  Cumberland  cultivar  with 
SOYGRO.  Thus,  model  traits  for  Williams  [MG  III]  were  used  initially  and 
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Table  6-1 .  Yield  data  for  Cumberland  cultivar--1 986  (Egli.1 988). 

Plants  m"2  Yield  (g  m"2)  Seeds  m*2  mg  seed'"' 

211 
205 
203 
204 
196 
206 
199 
200 
203 
191 
195 
197 
189 
192 
194 
195 
197 


0.63 

94.2 

442 

0.75 

136.4 

665 

0.90 

157.5 

777 

1.07 

198.5 

973 

1.27 

200.4 

1021 

1.51 

275.0 

1317 

1.08 

289.7 

1449 

2.14 

314.5 

1576 

2.56 

353.4 

1763 

3.04 

327.0 

1710 

3.62 

399.0 

2034 

4.32 

408.4 

2072 

5.14 

386.5 

2049 

6.12 

409.5 

2131 

7.29 

457.5 

2349 

8.69 

455.1 

2338 

12.32 

468.4 

2390 

Table  6-2.  Irrigation  schedule  for  Cumberland  cultivar-1 986  (Egli,1 988). 


Date  Rate  (mm) 


June 

17 

6.4 

June 

23 

12.7 

July 

28 

38.1 

August 

18 

25.4 
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Table  6-3.  Definitions  of  soil  parameters  and  characteristics  of  Maury  silt  loam. 


SALB- 

DMOD. 

DLAYR 

5 

15 

13 

12 

13 

23 

21 

17 

18 

26 

28 

40 

28 


0.13 
=  1.0 


U  -  9.78 
SWCON1 


0.00132 


SWCON  -  0.60 
SWCON2  -  82.3 


CN2  >  84.0 
SWCON3  -  6.67 


LL 

0.550 

0.550 

0.168 

0.184 

0.184 

0.203 

0.215 

0.227 

0.240 

0.248 

0.240 

0.279 

0.241 


0.200 
0.200 
0.308 
0.322 
0.322 
0.339 
0.349 
0.359 
0.368 
0.374 
0.368 
0.398 
0.368 


SAT 

0.390 

0.390 

0.391 

0.391 

0.391 

0.394 

0.395 

0.396 

0.403 

0.406 

0.404 

0.413 

0.406 


SW 

0.200 

0.200 

0.308 

0.322 

0.322 

0.339 

0.349 

0.359 

0.368 

0.374 

0.368 

0.398 

0.368 


1.000 
1.000 
1.000 
0.750 
0.750 
0.350 
0.200 
0.200 
0.200 
0.150 
0.100 
0.000 
0.000 


BQ 

1.43 
1.43 
1.43 
1.43 
1.43 
1.42 
1.42 
1.41 
1.39 
1.39 
1.39 
1.50 
1.39 


TAV  -  20.0 
RWUMX  -  0.04 

D£ 

2.14 

2.14 

0.99 

0.48 

0.48 

0.38 

0.36 

0.34 

0.38 

0.34 

0.34 

0.34 

0.34 


AMP  -  0.0 
PHAC3-1.0 

EH 

6.7 

6.7 

7.0 

6.4 

6.4 

6.0 

5.8 

5.5 

5.1 

5.5 

5.4 

5.1 

4.7 


SALB  Bare  soil  albedo  (no  units). 

U  Upper  limit  of  stage  1  soil  evaporation  (mm). 

SWCON        Soil  water  drainage  constant,  fraction  drained  per  day. 

CN2  SCS  curve  number  used  to  calculate  daily  runoff. 

TAV  Annual  average  ambient  temperature  ( °  C). 

At^P  Annual  amplitude  in  mean  monthly  temperature  ( °  C). 

DMOD  Zero-to-unity  factor  w/hich  reduces  the  rate  constant  for  mineralization  of  the  humus 

pool  for  soils  which  are  poor  mineralizers  due  to  chemical  or  physical  protection  of  the 

organic  matter  (default  =  1). 
SWC0N1       Coefficient  in  the  steady  state  solution  to  the  radial  flow,  root  uptake  equation,  cm  ^  / 

cm  root-day  (default  =  0.00132). 
SWCON2      CoeHicient  in  the  steady  state  solution  to  the  radial  flow,  root  uptake  equation,  cm  ^  / 

cm  root -day  (value  calculated  by  retrieval  program). 
SWC0N3      Coefficient  in  the  steady  state  solution  to  the  radial  flow,  root  uptake  equation,  cm  ^  / 

cm  root-day  (default  =  6.67). 
RWUMX        Maximum  daily  root  water  uptake  per  unit  root  length,  cm  ^/  cm  root-day  (default  = 

0.04). 
PHFAC3        Variable  to  reduce  apparent  photosynthesis  attributed  to  soil  fertility  (for  grain  legume 

nx)dels,  default  =  1.00). 
DLAYR  (L)     Thickness  of  soil  layer  L  ,  cm. 

LL  (L)  Lower  limit  of  plant-extractable  soil  water  for  soil  layer  L,  cnr/crrr. 

DUL  (L)  Drained  upper  limit  soil  water  content  for  soil  layer  L,  crrr/civr. 

SAT  (L)  Saturated  water  content  for  soil  layer  L,  cm^/cm^ 

SW  (L)  Default  soil  water  content  for  soil  layer  L,  cnv^/crrr. 

WR  (L)  Weighting  factor  for  soil  depth  L  to  determine  new  root  growth  distribution,  no  units. 

BD  (L)  Moist  bulk  density  of  soil  in  soil  layer  L,  g/crr?. 

OC  (L)  Organic  cartx)n  concentration  in  soil  layer  L,  %. 

PH  (L)  pH  of  soil  in  soil  layer  L  in  a  1 :1  soil:  water  slurry. 
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modified.     The  genetic  and  plant  parameters  which  were  altered  in  model 
calibration  were: 


SHVAR: 


SDVAR: 


Maximum  rate  of  dry  matter  accumulation  per  shell  (pod  walls)  at 
optimum  temperature  and  assimilation  supply. 

Maximum  rate  of  dry  matter  accumulation  per  seed  during  the 
linear  filling  period  and  under  optimum  temperature  and 
assimilate  supply. 


VARTH  (6): 

VARTH(10): 

PODVAR: 

FLWVAR: 

SCV: 

LFMAX: 

YVSHT: 

YVSWH: 


Photothermal  time  from  R1  to  the  time  of  beginning  of  podset. 

Photothermal  time  from  R1  to  R7  (beginning  maturity). 

Maximum  number  of  pods  produced  per  m^  per  day  under  short 
days  and  optimum  temperature  and  assimilate  supply. 

Maximum  number  of  flowers  produced  per  m^  per  day  under 
short  days  and  optimum  temperature  and  assimilate  supply. 

Scattering  coefficient  for  light  (fraction  reflected  and  transmitted 
by  a  single  leaf). 

Light-saturated  photosynthetic  rate  of  an  upper  mature  leaf  at 
optimum  temperature,  N  and  water  status. 

Parameters  which  determine  canopy  height  via  expected 
internode  length  as  a  function  of  vegetative  stage. 

Paramaters  which  determine  canopy  width  via  petiole  and  leaf 
extension  as  a  function  of  vegetative  stage. 


Beginning  with  genetic  and  crop  traits  of  Williams  cultivar  [MG  III],  the  R1 
date  was  predicted  within  1  day,  but  seed  size  and  seed  yield  were  too  low. 
Since  Egli  had  not  observed  the  R7  date  and  because  SOYGRO  had  previously 
predicted  too  early  maturity  for  the  Williams  cultivar  at  the  Lexington  site,  genetic 
traits  (VARTH  (10))  were  altered  to  increase  the  seed  fill  period  ( time  to  R7  )  and 
increase  seed  size  (SHVAR  and  SDVAR).    Varth  (6)  was  reduced  to  initiate 
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podset  sooner  after  R1.  LFMAX,  PODVAR  and  FLWVAR  were  all  increased  to 
increase  seed  yield  (Table  6-4.  ).  Increasing  LFMAX  has  the  effect  of  increasing 
seed  yield  about  4.6%  per  10%  increase  in  LFMAX.  The  SCV  value  was 
reduced  to  increase  light  capture,  especially  at  higher  plant  densities,  in  order  to 
increase  light  interception  to  values  comparable  to  limited  data  from  Egli  (1988). 

Length  per  internode  was  reduced  from  6.5  cm/internode  to  5.3 
cm/internode.  Petiole-leaflet  extension  (plant  diameter)  per  vegetative  stage  was 
reduced  from  6.7  to  5.3  cm.  The  reduction  of  the  height  and  width  of  the  canopy 
hedgerow  was  done  to  reduce  light  absorption  and  hence  photosynthesis  of 
widely  spaced  plants  relative  to  higher  plant  densities  (closed  canopies).  This 
calibration  resulted  in  a  plant  simulated  height  and  width  of  80  and  70  cm, 
respectively  at  maximum  size  and  a  vegetative  stage  of  16.7.  This  calibration 
influenced  the  shape  of  the  yield  response  to  density,  primarily  the  rapid  yield 
increase  phase. 

Results  of  Model  Calibration 

Model  calibration  was  conducted  using  Egli's  (1988)  yield  data. 
Simulated  yield  data  is  provided  in  Table  6-5.  Simulated  seed  weight  (g  m'^) 
data  generated  by  SOYGRO  approximated  to  Egli's  (1988)  observed  data  (Fig. 
6-1).  However,  there  was  a  tendency  for  SOYGRO  to  underestimate  seed 
weight  data  at  plant  densities  greater  than  two  plants  m"^.  Simulated  and 
observed  weight  seed  "''  showed  opposite  trends;  the  former  increased  while  the 
latter  decreased  with  plant  density  (Fig.  6-2).  However,  average  weights  seed'^ 
for  simulated  data  and  observed  data  approximated  to  each  other.  Seed  number 
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Table  6-4.  Summary  of  changes  in  genetic  parameters  of  Williams 
cultivar  and  cropparameters  in  calibrating  SOYGRO  using  Cumberland 
cultivar,1986. 


Genetic  parameters 

Williams  "initial" 

Cumberland 

cultivar 

cultivar 

SHVAR  (rngsheir""  d"""): 

11.0 

13.4 

SDVAR  (mg  seed"""  d"""): 

6.0 

6.9 

VARTH(6)(PTdays*): 

6.0 

4.0 

VARTH(10)(PTdays): 

44 

48 

PODVAR  (pods  m-2  d"""  ): 

120 

165 

FLWVAR  (flowers  m-2  d"""  ): 

240 

330 

LFMAX(mgC02m-2s-''): 

1.046 

1.287 

Crop  parameters 

Initial 

New 

SCV  (fraction): 

0.20 

0.15 

YVSHT(cm  internode"''): 

6.5 

5.3 

YVSWH  (cminternode"''): 

6.7 

5.3 

*  photothermal  days 
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Table  6-5.  Simulated  yield  data  for  Cumberland  cultivar,1 986 


Plants  m"2 

Yield  (g  m'^) 

Seeds  m"2 

mg  seed""" 

0.63 

113.1 

654 

173 

0.75 

133.3 

769 

173 

0.90 

157.7 

907 

174 

1.07 

184.2 

1056 

175 

1.27 

213.5 

1215 

176 

1.51 

246.0 

1368 

180 

1.08 

280.0 

1565 

179 

2.14 

303.8 

1622 

190 

2.56 

334.3 

1714 

195 

3.04 

355.8 

1796 

198 

3.62 

375.2 

1882 

199 

4.32 

392.1 

1946 

201 

5.14 

407.0 

2015 

202 

6.12 

418.8 

2049 

204 

7.29 

429.3 

2085 

206 

8.69 

442.2 

2156 

205 

12.32 

462.8 

2268 

204 
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Fig.  6-1.  Observed  and  simulated  effects  or  plant  population 
density  on  seed  yield  of  Cumberland  cultivar  at  maturity^  986. 
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Fig.  6-2.  Observed  and  simulated  effecbs  of  plant  population  density 
0  n  we ig ht  per  seed  of  C umberia  nd  c  u  Itivar  at  maturity,  1986. 
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(m"2)  data  behaved  in  a  manner  similar  to  seed  weight  (g  m'2),  but  there  was  a 
tendency  for  SOYGRO  to  underestimate  seed  number  (m'^)  data  at  plant 
densities  greater  than  two  plants  m"^  (Fig.  6-3).  Simulated  yield  per  plant  data 
closely  approximated  observed  data  (Fig.  6-4). 

An  Investigation  of  Duncan's  Assumptions 

Duncan  considered  percent  light  interception  at  R5  as  an  estimate  of  the 
amount  of  light  intercepted  by  the  crop  between  R5  and  R7  (physiological 
maturity).  However,  Fig.  6-5.  shows  that  percent  light  intercepted  at  R5  (  noon 
and  daily)  are  not  true  measures  of  light  intercepted  between  R5  and  R7;  it  is  only 
cumulative  percent  light  intercepted  between  R5  and  R7  which  is  a  direct 
measure  of  light  intercepted  during  the  reproductive  period.  The  greater  percent 
light  interception  observed  for  the  R5  to  R7  period  may  be  due  to  the  fact  that  for 
indeterminate  cultivars  SOYGRO  predicts  R5  from  the  beginning  of  seed 
initiation  at  lower  nodes;  however,  beginning  of  seed  initiation  occurs  prior  to  the 
end  of  leaf  expansion,  hence  a  greater  percent  of  light  will  be  intercepted  by  the 
crop  after  the  beginning  of  seed  growth  on  indeterminates.  The  problem  stems 
from  the  R5  definition  of  Fehr  et  al  (1971).  Although  R5  is  called  "beginning 
seed",  in  fact  the  stage  is  observed  only  at  the  upper  four  expanding  nodes,  an 
dactual  seed  growth  was  initiated  several  weeks  earlier  at  lower  nodes  on 
indeterminate  plants.  For  determinate  cultivars,  this  problem  is  less  acute. 

From  Duncan's  perspective,  remobilization  contributes  to  final  yield. 
However,  in  estimating  changes  in  seed  yield  efficiency  with  plant  density, 
remobilization  was  disregarded  as  a  process  which  can  account  for  changes  in 
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Fig.  6-3.  Observed  and  simulated  effects  of  plant  population  density 
on  seed  nunnberper  unitareaof  Cunnberlandcultivarat  maturity,  1986. 
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Fig.  6-4.  Obsep/ed  and  simulated  effects  of  plant  population 
density  on  seed  yield  per  plant  of  Cumberland  cultivar  at 
maturity,  1 986. 
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Fig.  6-5.  The  relationships  between  light  intercepted  betv/een  R5 
and  R7,  and  percent  light  absorbed  at  noon  (R5),  daily  (R5)  and 
between  R5  and  R7  for  Cumberland,  1 986. 
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the  efficiency  with  which  inputs  are  used  to  produce  seed.  In  Fig.  6-6.,  the  total  N 
per  unit  area  mobilized  increases  with  plant  population  density.  This  result  is 
expected  because  at  higher  plant  population  densities  there  are  more  nutrients 
per  unit  area  for  remobilization.  However,  the  extent  of  remobilization,  as 
estimated  by  the  percent  of  mobilized  nitrogen  in  seeds,  also  tends  to  increase 
with  plant  density  (Fig.  6-7.).  To  the  extent  that  the  percentage  of  mobilized  N  in 
seed  is  increased,  that  fraction  of  seed  was  produced  with  less  photosynthate 
since  the  cost  of  N2  fixation  was  paid  for  during  vegetative  growth  and  some  C 
atoms  were  remobilized  with  the  amino  acids.  The  effect  would  be  increased 
"apparent"  light  use  efficiency.  Therefore,  the  process  of  remobilization  will 
confound  comparisons  of  seed  yield  efficiencies  made  between  crops  of  different 
plant  densities.  If  percent  remobilized  N  in  the  seed  remained  constant  with 
changes  with  plant  density,  the  remobilization  would  impact  only  upon  yield 
differences  and  not  differences  in  seed  yield  efficiency. 

Based  on  the  results  of  SOYGRO,  Duncan's  first  two  assumptions  do  not 
hold  and  hence  the  following  are  clear.  First,  it  is  inappropriate  to  use  the  light 
interception  at  R5  to  calculate  the  efficiency  with  which  light  is  used  to  produce 
seed  between  R5  and  R7.  Secondly,  remobilization  has  a  differential  effect  on 
seed  yield  and  therefore  confounds  calculations  of  seed  yield  efficiency. 
Therefore,  seed  yield  efficiencies  using  light  intercepted  between  R5  and  R7  are 
invalid  unless  the  differential-remobilization-effect  is  removed.  In  calculating  and 
comparing  the  efficiency  with  which  light  is  used  by  crops  grown  at  different  plant 
densities  three  approaches  will  be  used:  total  light  use  efficiency-glucose 
equivalent,  seed  yield  efficiency  and  seed  yield  efficiency-glucose  equivalent. 
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Fig.  6-6.  The  simulated  effect  of  plant  population  density  on  N 
mob  ilizedperun  it  area  forCumberiandcultivar, 1986. 
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Fig.  6-7.  The  simulated  effect  of  plant  population  density  on 
percent  mob i I ized-N  in  seeds  of  Cumberland  cultivar,  1986. 
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Total  Light  Use  Efficiencv--Glucose  Equivalent 

Total  light  use  efficiency  of  a  crop  is  determined  by  converting  the  total 
biomass  produced  (including  abscised  organs)  by  a  crop  into  its  glucose 
equivalent,  and  dividing  this  value  by  the  total  light  intercepted  by  the  crop.  It  is 
assumed  that  the  weight-ratio  of  vegetative  parts  of  the  soybean  crop  is  not 
affected  by  plant  density.  This  ratio  for  stem  to  leaf  to  shell  to  root  which  is 
derived  from  SOYGRO  is  1.87  to  1.15  to  1.02  to  1.00,  respectively  (Jones  et 
al.,1989).  It  is  also  assumed  that  to  produce  one  gram  of  stem,  leaf,  shell  and 
root  tissue,  1.44,  1.61,  1.49  and  1.37  grams  of  glucose  are  required,  respectively. 
Therefore,  the  average  cost  in  glucose  per  gram  of  vegetative  tissue  produced  is 
1.49  g.  Seed  production  is  assumed  to  cost  2.08  g  of  glucose  for  each  gram  of 
seed  produced  with  fixed  N,  and  1.45  g  of  glucose  for  each  gram  of  seed 
produced  with  mobilized  N.  The  effect  of  plant  density  on  total  and  seed  yield 
light  use  efficiency  is  given  in  Table  6-6  and  Fig.  6-8.  Between  populations  of 
0.63  and  1.08  plants  m  '^,  where  interplant  competition  is  negligible,  it  appeared 
that  total  light  use  efficiency  was  constant.  The  increase  in  total  light  use 
efficiency  with  increasing  plant  density  is  caused  by  a  greater  frcation  of  the  leaf 
area  operating  at  lower  light  intensities  where  their  quantum  efficiency  is  greater. 

Seed  Yield  Efficiency  and  Seed  Yield  Efficiencv-Glucose  Equivalent 

Both  seed  yield  efficiency  and  seed  yield  efficiency-glucose  equivalent 
were  determined  by  allowing  for  the  mobilization  of  nutrients  from  vegetative 
tissue.  However,  if  a  "true"  seed  yield  efficiency  can  be  computed,  these  nutrients 
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Table  6-6.  Total  light  use  efficiencies  as  determined  by  SOYGRO  for 
Cumberland  cultivar,1986  at  seventeen  plant  densities. 


Glucose  equivalents  (g  glucose  m  "^) 

Total  light 

Plant 

Vegetation 

Seed 

Total 

Total  light  use 

density 

absorbed 

efficiency 

Plants  m"2 

(MJ  m  -2) 

(g  glucose  MJ 

0.63 

204.9 

214.5 

419.4 

601 

0.698 

0.75 

241.8 

252.6 

494.3 

707 

0.699 

0.90 

285.8 

298.2 

584.0 

836 

0.699 

1.07 

334.0 

348.4 

682.5 

976 

0.699 

1.27 

338.6 

404.5 

793.1 

1131 

0.701 

1.51 

448.8 

465.3 

914.1 

1303 

0.702 

1.08 

512.7 

529.9 

1042.5 

1486 

0.702 

2.14 

597.7 

573.3 

1171.0 

1653 

0.708 

2.56 

663.8 

629.1 

1292.9 

1806 

0.716 

3.04 

725.5 

669.0 

1394.4 

1933 

0.721 

3.62 

776.9 

703.3 

1480.2 

2047 

0.723 

4.32 

826.1 

734.5 

1560.6 

2150 

0.726 

5.14 

873.8 

760.4 

1634.2 

2244 

0.728 

6.12 

923.1 

780.2 

1703.3 

2332 

0.730 

7.29 

968.4 

797.8 

1766.3 

2413 

0.732 

8.69 

999.8 

821.0 

1820.8 

2491 

0.731 

12.32 

1073.2 

854.9 

1928.0 

2640 

0.730 
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Fig.  6-8.  The  simulated  effect  or  plant  population  density  on 
total  light  use  efficiency  of  Cumberland  cultivar,  1 986. 
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Table  6-7.  Seed  yield  efficiencies  as  determined  by  SOYGRO  for  Cumberland 
cultivar,1986  at  seventeen  plant  densities. 


No  mobilization 

Mobilization 

Plant 

Seed  dry 

Light 

Seed  yield 

Seed  yield 

Seed  yield 

density 

weight 

absorbed 

efficiency 

efficiency 

efficiency 

Plants  m'2 

(gm-2) 

(R5  -  R7} 

(gMJ-l) 

(g  glucose  MJ  ''^) 

(gMJ-1) 

(MJm 

-^) 

0.63 

106.0 

515 

0.206 

0.428 

0.219 

0.75 

124.8 

605 

0.206 

0.429 

0.220 

0.90 

147.5 

714 

0.207 

0.430 

0.220 

1.07 

171.9 

831 

0.207 

0.430 

0.221 

1.27 

198.2 

960 

0.207 

0.430 

0.222 

1.51 

227.4 

1101 

0.207 

0.430 

0.223 

1.08 

257.2 

1247 

0.206 

0.429 

0.224 

2.14 

284.4 

1373 

0.207 

0.431 

0.224 

2.56 

308.2 

1477 

0.209 

0.434 

0.226 

3.04 

327.9 

1554 

0.211 

0.439 

0.229 

3.62 

346.0 

1616 

0.214 

0.445 

0.232 

4.32 

360.2 

1667 

0.216 

0.449 

0.235 

5.14 

374.4 

1710 

0.219 

0.456 

0.238 

6.12 

385.4 

1745 

0.220 

0.459 

0.239 

7.29 

396.0 

1778 

0.223 

0.463 

0.241 

8.69 

404.5 

1805 

0.224 

0.466 

0.245 

12.32 

420.4 

1852 

0.227 

0.472 

0.251 
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Fig.  6-9.  The  simulated  effects  of  plant  population  density 
on  seed  yield  efficiency  of  Cumberland  cultivar,  1 986,  with 
and  without  the  benefit  of  protein  remobillzatlon. 
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must  be  theoretically  discarded  and  not  put  into  seed  production.  The  resulting 
seed  yield  achieved  (strictly  the  result  of  photosynthetic  activity  between  R5  and 
R7)  was  then  divided  by  the  light  intercepted  by  the  crop  between  growth  stages 
R5  and  R7.  The  results  as  hypothesized,  are  provided  in  Table  6-7.  and  Fig.  6.9. 
With  no  protein  mobilization,  "true"  seed  yield  efficiency  increased  with  plant 
population  density  for  both  harvested  seed  material  and  glucose  equivalents. 
Apparent  seed  yield  efficiency  was  increased  as  a  result  of  protein  remobilization 
at  all  plant  densities,  but  the  effect  was  greater  as  plant  density  increased. 

Conclusions 

Results  of  SOYGRO  showed:  i)  Duncan's  assumptions  that  light 
interception  at  R5  is  an  effective  basis  for  calculating  and  comparing  the  seed 
yield  efficiencies  of  crops  grown  at  different  plant  population  densities  does  not 
appear  to  hold;  ii)  remobilization  is  affected  by  plant  density  and  therefore,  has  a 
differential  impact  on  seed  yield.  As  a  result,  calculations  of  seed  yield  efficiency 
to  accurately  estimate  a  crop's  productive  capacity  after  R5  (beginning  seed) 
must  remove  the  effect  of  remobilization  on  seed  yield. 

SOYGRO  estimations  of  total  light  use  efficiency  (glucose  equivalents) 
and  seed  yield  efficiency  (remobilization  effect  removed),  indicated  that  Duncan's 
assumption  that  with  increasing  plant  population  there  is  an  increased  efficiency 
with  which  light  is  used  for  seed  production,  appears  to  hold. 

SOYGRO  holds  much  potential  for  further  investigations  of  yield-plant 
density  responses.   Future  studies  should  focus  on  the  morphological  or 
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physiological  basis  for  the  increased  total  and  seed  yield  efficiencies  predicted 
by  SOYGRO. 


CHAPTER  7 
SUMMARY  AND  CONCLUSIONS 


Experimental  results  show  that  density-mediated  plastic  responses  and 
genotype  play  a  major  role  in  the  process  of  yield  achievement  in  soybean 
[Glycine  max  (I.)  Merr.].  Preliminary  investigations  into  yield-plant  density 
responses  with  SOYGRO-a  process  oriented  crop  growth  simulator  for  soybean, 
show  that  one  consequence  of  density-mediated  plastic  responses  is  an 
increased  efficiency  with  which  intercepted  light  is  used  for  seed  production. 

In  the  1987  field  experiment,  plant  population  had  no  significant  effect  on 
the  variables  measured  at  R5,  R6  and  R8.  However,  the  effect  of  cultivar  on 
soybean  traits  allowed  for  yield  investigations.  Correlation  analyses  on  the 
data  sets  for  each  cultivar,  and  for  the  combined  data  set  of  all  cultivars,  at  both 
R6  and  R8,  showed  that  pod  dry  weight  (  g  m  '^),  seed  number  (m"2),  and 
pod  number  (m  '^)  were  all  highly  correlated  with  seed  dry  weight  (g  m"^).  Seed 
number  pod'"'  was  unaffected  by  plant  population  density  and  Kirby  had  the 
lowest  value.  However  the  greater  weight  seed  ""•  (g)  of  Kirby  contributed  to 
Kirby    having  a  greater    seed    dry  weight     (g  m  '^)  at  R8  than  Forrest. 

Weight  seed  "^  (g)  was  neither  plastic  nor  was  it  highly  correlated  with 
seed  dry  weight  (g  m"^).  However,  among  yield  components,  weight  seed  '^ 
accounted  for  the  yield  difference  observed  among  cultivars.   Being  under 
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genetic  control,  weight   seed  "''  can  account  for  differences  In  genotypes,  but  is 
unable  to  account  for  density-mediated  plastic  response  of  seed  yield. 

In  the  1988  field  experiment  ,  the  variables  measured  at  growth  stages 
R3,  R5  and  R8  included  total  dry  weight,  leaf,  stem  and  pod  fractions,  leaf  area 
index,  light  interception,  N  and  P  concentrations  of  leaf  tissue,  pod  dry  weight  and 
number,  and  seed  dry  weight  and  number.  It  was  shown  that  all  reproductive, 
photosynthetic  and  remobilization  parameters  which  were  plastic  were 
associated  with  seed  yield.  The  variables  most  sensitive  to  plant  population,  and 
which  were  also  closely  associated  with  yield  were  seed  number  and  pod 
number  at  R8  and  total  dry  weight  and  leaf  dry  weight  at  R5.  Nutrient  content 
(g  m'2)  was  moderately  associated  with  yield;  however,  nutrient  concentration 
was  neither  associated  with  yield  nor  was  it  affected  by  plant  density. 

Results  of  SOYGRO  showed:  i)  Duncan's  assumptions  that  light 
interception  at  R5  is  an  effective  basis  for  calculating  and  comparing  the  seed 
yield  efficiencies  of  crops  grown  at  different  plant  population  densities  does  not 
appear  to  hold;  ii)  remobilization  is  affected  by  plant  density  and  therefore,  has  a 
differential  impact  on  seed  yield.  As  a  result,  calculations  of  seed  yield  efficiency 
which  estimate  a  crop's  productive  capacity  after  R5  (beginning  seed)  must 
remove  the  effect  of  remobilization  on  seed  yield  to  be  accurate  . 

SOYGRO  estimations  of  total  light  use  efficiency  (glucose  equivalents) 
and  seed  yield  efficiency  (remobilization  effect  removed),  indicated  that  Duncan's 
assumption  that  with  increasing  plant  population  there  is  an  increased  efficiency 
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with  wNch  light  is  used  for  seed  production,  appears  to  hold.  The  remobilization 
effect  adds  to  the  increase  in  apparent  seed  yield  efficiency  with  increasing  plant 
population. 

Further  investigations  into  the  ecological  and  physiological  mechanisms 
by  which  plant  population  density  affects  seed  yield,  appeared  to  require  the 
collection  of  data  over  the  entire  growing  season.  Such  data  should  include 
plant  demographic  variables  such  as  leaf,  pod,  branch  number  and  distribution; 
also  needed  are  seasonal  light  interception  data.  SOYGRO  has  an  important 
role  to  play  in  developing  hypotheses  and  guiding  research  efforts  into 
yield-plant  density  responses. 
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